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Abstract. In the paper we review some recent results of the theory of hierarchies 
of quantum evolution equations. 

The evolution of infinite-particle quantum systems is described within the framework 
of the evolution of states by the quantum BBGKY hierarchy for marginal density 
, operators or within the framework of an equivalent approach in terms of the marginal 

Q^l observables by the dual quantum BBGKY hierarchy. The nonperturbative solutions 

of the Cauchy problem of the dual quantum BBGKY hierarchy and the quantum 
BBGKY hierarchy are constructed. 

The origin of the microscopic description of non-equilibrium quantum correlations 
is considered. It consists in one more approach to the description of the evolution 
of states of quantum many-particle systems by means of correlation operators. The 
correlation operators and the marginal correlation operators are governed by the 
' von Neumann hierarchy and the quantum nonlinear BBGKY hierarchy respectively. 

. The nonperturbative solutions of the Cauchy problem of both these hierarchies are 

constructed. 

We develop also an approach to a description of the evolution of states by means 
of the quantum kinetic equations. For initial states which are specified in terms of 
a one-particle density operator the equivalence of the description of the evolution 
of quantum many-particle states by the Cauchy problem of the quantum BBGKY 
hierarchy and by the Cauchy problem of the generalized quantum kinetic equation 
^ . together with a sequence of explicitly defined functionals of a solution of stated 

^ I kinetic equation is established. The relationships of the specific quantum kinetic 

equations with the generalized quantum kinetic equation are discussed. 
In conclusion the mean field asymptotic behavior of stated hierarchy solutions are 
established. The constructed asymptotics are governed by the quantum Vlasov hier- 
archy for limit states, the nonlinear quantum Vlasov hierarchy for limit correlations 
and the dual quantum Vlasov hierarchy for limit observables respectively. 
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1 Introduction: motivations and results 

Experimental advances in the Bose condensation of atomic gases and in the strong correlated 
Fermi systems have stimulated interesting problems in the theory of evolution equations of 
quantum many-particle systems. Among them it is a description of collective behavior of 
interacting particles by quantum kinetic equations, i.e. the evolution equations for a one- 
particle marginal density operator |ll|2]. The rigorous derivation of quantum kinetic equations 
is a very trendy subject nowadays, and the number of publications devoted to it has grown 
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last decade. In particular the considerable progress in the rigorous derivation of the nonlinear 
Schrodinger equation and the Gross- P it aevskii equation for the Bose condensate [8HT9] in mean 
field scaling limit as well as the quantum Boltzmann equation [20l[2T] is observed. Owing to 
the intrinsic complexity and richness of this problem, first of all it is necessary to develop 
an adequate mathematical theory of evolution equations of quantum many-particle systems 
underlying of kinetic equations. 

In this review we expound the foundation of the theory of quantum evolution equations of 
statistical mechanics. We construct nonperturbative solutions of the Cauchy problem of the 
quantum BBGKY hierarchy for marginal density operators and the dual quantum BBGKY 
hierarchy for marginal observables on appropriate Banach spaces. On basis of obtained results 
the origin of the microscopic description of non-equilibrium quantum correlations is considered. 
In particular, the nonperturbative solutions of the Cauchy problem of the von Neumann hi- 
erarchy for correlation operators and the nonlinear quantum BBGKY hierarchy for marginal 
correlation operators are constructed. We develop also one more approach to the description of 
the evolution of quantum many-particle systems by means of the quantum kinetic equations. 
For initial states which are specified in terms of a one-particle density operator the equivalence 
of the description of the evolution of quantum infinite-particle states by the Cauchy problem 
of the quantum BBGKY hierarchy and by the Cauchy problem of the generalized quantum 
kinetic equation is established. We also describe the mean field scaling asymptotic behavior of 
nonperturbative solutions of hierarchies of quantum evolution equations under consideration. 

We outline the structure of the paper and the main results. In introductory section 2 we 
set forth the traditional approach to the description of the evolution of quantum many-particle 
systems. It is well known that a description of quantum many-particle systems is formulated 
in terms of two sets of objects: observables and states. The functional of the mean value of 
observables defines a duality between observables and states and as a consequence there exist 
two equivalent approaches to the description of the evolution, namely in terms of the evolution 
equations for observables (the Heisenberg equation) and for states (the von Neumann equation). 
We adduce some preliminary facts about dynamics of finitely many quantum particles described 
by these evolution equations within the framework of nonequilibrium grand canonical ensemble. 

In section 3 one more an equivalent approach to the description of the evolution of states 
of quantum many-particle systems is given by means of operators which are interpreted as 
correlation operators. To justify the von Neumann hierarchy for correlation operators, at 
first we consider in detail the motivation of the introduction of correlation operators to the 
description of states or in other words, the origin of the microscopic description of correlations 
in quantum many-particle systems is considered. 

In section 4 for the description of the evolution of infinite-particle quantum systems we 
introduce the hierarchies of evolution equations for the marginal observables, density operators 
and correlation operators. They are derived as the evolution equations for one more but an 
equivalent method of the description of states and observables of finitely many particles. The 
nonperturbative solutions of the Cauchy problem of the dual quantum BBGKY hierarchy for 
the marginal observables, the quantum BBGKY hierarchy for marginal density operators and 
the nonlinear quantum BBGKY hierarchy for marginal correlation operators are constructed 
as an expansion over particle clusters which evolution is governed by the corresponding-order 
cumulant (semi-invariant) of the groups of operators of the Heisenberg equations, the von 
Neumann equations and the von Neumann hierarchy of finitely many particles respectively. 
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In section 5 we develop an approach to a description of the evolution of states by means of 
the quantum kinetic equations. For initial states which are specified in terms of a one-particle 
density operator the equivalence of the description of the evolution of quantum many-particle 
states by the Cauchy problem of the quantum BBGKY hierarchy and by the Cauchy problem 
of the generalized quantum kinetic equation together with a sequence of explicitly defined 
correlation functionals of a solution of stated quantum kinetic equation is established. The 
relationship of the generalized quantum kinetic equation with the specific quantum kinetic 
equations is discussed. 

In section 6 the mean field (self-consistent field) asymptotic behavior of constructed above so- 
lutions of the quantum hierarchies and the generalized quantum kinetic equation is established. 
In particular, the quantum Vlasov hierarchy for limit states, the nonlinear quantum Vlasov 
hierarchy for limit correlations and the dual quantum Vlasov hierarchy for limit observables 
are rigorously derived. 

Finally in section 7 we conclude with some observations and perspectives for future research. 

2 On evolution equations of quantum many-particle sys- 
tems 

It is well known that a description of quantum many-particle systems is formulated in terms 
of two sets of objects: observables and states. The functional of the mean value of observables 
defines a duality between observables and states and as a consequence there exist two approaches 
to the description of the evolution, namely in terms of the evolution equations for observables 
(the Heisenberg equation) and for states (the von Neumann equation). In this section we 
adduce some preliminary facts about dynamics of finitely many quantum particles described 
within the framework of nonequilibrium grand canonical ensemble. 

2.1 Preliminaries 

We consider a quantum system of a non-fixed, i.e. arbitrary but finite, number of identical 
(spinless) particles (nonequilibrium grand canonical ensemble) obeying the Maxwell-Boltzmann 
statistics in the space W. We will use units where h = 21x1% = 1 is a Planck constant, and 
m = 1 is the mass of particles. 

Let "H be a one-particle Hilbert space, then the n-particle space "Hn, is the tensor product of n 
Hilbert spaces "H. We adopt the usual convention that l-i^^ = C. We denote by J-'-u = 0^o'H„ 
the Fock space over the Hilbert space "H. 

Let a sequence g = {go,gi, . . . ,gn, ■ ■ ■) be an infinite sequence of self-adjoint bounded op- 
erators defined on the Fock space J-'-h and go G C. The operator gn defined on the n-particle 
Hilbert space Tin = 'H®'^ will be denoted by gn{^, ■ ■ ■ ,n). Let the space -C(J-'-h) be the space of 
sequences g = {go, gi, . . . , gn, ■ ■ ■) of bounded operators gn defined on the Hilbert space Tin that 
satisfy symmetry condition: (^^(l, . . . ,n) = gnih, • • • , ^n), for arbitrary (ii, . . . , z„) G (1, . . . , n), 
equipped with the operator norm ||.||£c^^) |5j. We will also consider a more general space 
ii^(J^'^) with a norm 
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where < 7 < 1. We denote by £^^o(-^-h) C 2,^{Fh) the everywhere dense set of finite 
sequences of degenerate operators with infinitely differentiable kernels with compact supports. 
Observables of finitely many quantum particles are sequences of self-adjoint operators from the 
space [5]. 

Let 2}{J^y) = ^'^=0^.^ i'Hn) be the space of sequences / = (/o, /i, ...,/„,. . .) of trace class 
operators /„ = /„(!, . . . ,n) G il^CHn) and /o G C, satisfying the mentioned above symmetry 
condition, equipped with the trace norm 

00 00 
IULi(-^h) ^ ll'^"LiCHn) ^ X] Tri^..._„|/„(1, . . . ,n)|, 

n=0 n=0 

where Tri^ „ is the partial trace over 1, . . . ,n particles. The everywhere dense set of finite 
sequences of degenerate operators with infinitely differentiable kernels with compact supports 
from the space £^{J-'-u) we denote by The sequences of self-adjoint operators G 

ii^(Hn), n > 1, which kernels are known as density matrices, defined on the n-particle Hilbert 
space Tin = i^^(M^"), describe states of a quantum system of non- fixed number of particles. 
The space 2:^{J^n) is dual to the space with respect to the bilinear form 

00 

{9, /) = -T Tri,...,„ gnfn, (1) 

n=0 

where gn G and /„ G The range of positive continuous linear functional ([1]) on 

the space of observables ^{J^u) is interpreted as the mean values (averages) of observables. 

The Bose-Einstein and Fermi-Dirac statistics endow observables and states with additional 
symmetry properties |B] in comparison with the Maxwell-Boltzmann statistics. Let H he a. 
one-particle Hilbert space, then the n-particle spaces "H^, are correspondingly symmetric and 
antisymmetric tensor products of n Hilbert spaces "H that are associated with ri-particle systems 
of bosons and fermions [6] . We denote by J-"^ = 0^q'H^ the Bose and Fermi Fock spaces over 
the Hilbert space T-L respectively. 

The symmetrization operator and the anti-symmetrization operator S~ on 7^®" are 
defined by the formula 

^n=^Y.(±lPP-, (2) 

where the operator is a transposition operator of the permutation tt from the permutation 
group &n of the set (1, . . . , n) and |7r| denotes the number of transpositions in the permutation 
71. The operators are orthogonal projectors, i.e. {S^Y = Si^, ranges of which are corre- 
spondingly the symmetric tensor product and the antisymmetric tensor product T-L~ of n 
Hilbert spaces "H. 

2.2 The Heisenberg equation: the evolution of observables 

The Hamiltonian Hn of n-particle system is a self-adjoint operator with the domain V{Hn) C 

n n 

Hn = Y.m+ J2 Hn,^2), (3) 

i=l ii<i2 = l 



Theory of quantum evolution equations 



6 



where K{i) is the operator of a kinetic energy of the i particle and $(21,22) is the opera- 
tor of a two-body interaction potential. The operator K{i) acts on functions ipn that be- 
long to the subspace LliW"^) C V{Hn) C L'^iW'^) of infinitely different iable functions with 
compact supports according to the formula: K{i)'i/jn = —^Ag.ipn- Correspondingly we have: 
$(21, i2)4'n = ^{Qh , Qi2)'^n, and we assume that the function ^{qi^ , qi^) is symmetric with respect 
to permutations of its arguments, translation-invariant and bounded function. 

The evolution of observables A{t) = {Aq, Ai{t,l), . . . , An{t, 1, . . . , n), . . .) is described by the 
initial-value problem of a sequence of the Heisenberg equations [HIS] 

jA{t)=^A{t), (4) 
Ait)\t=o = AiO), (5) 

where A{0) = {Ao,Al{l), . . . , . . . . . .) G if ^ G V{Af) C the generator 

M = 0^qA/'„ is defined by the formula 

A/'n gn = -i{gnHn - Hngn) , (6) 

and Hn is the Hamiltonian ([3]). 

To determine a solution of the Cauchy problem dl])-® we introduce on the space ii(J^-H) the 
following one-parameter family of operators Q{t) = ©^o^^il^)- 

Qnit)gn = e''''-gne-''''-. (7) 

On the space the one-parameter mapping 3 t h-t- Q(t)g defines an isometric *- 

weak continuous group of operators, i.e. it is a Cg-group [5l|6]. The infinitesimal generator 
Af = '^^ ^his group of operators is a closed operator for the *-weak topology and on 

its domain of definition V{Afn) C QiTin) which is everywhere dense for the *-weak topology, 
Mn is defined as follows in the sense of the *-weak convergence of the space 

w*- lim ^(^„(t)^„ - gn) = -i{gnHn - Hngn), (8) 

where the operator Mngn = —iignHn — Hngn) is defined on the domain V[Hn) C T-Ln- 

The group of operators ^ preserves the self-adjointness of operators and the canonical 
commutation relations which are determined the physical meaning of operators characterizing 
particles and the structure of generator ([6]) of the Heisenberg equation (jl]). 

On the space for initial- value problem dl])-© the following statement holds. 

Theorem 1. A unique solution of initial-value problem dlD-Q of the Heisenberg equation is 
determined by the formula 

Ait) = g{t)A{0), (9) 

where the one-parameter family {^(t)}teiR of operators is defined by expression ([7j). For A{0) G 
V{Af) C Sy-yiTy), operator (Q is a classical solution and for arbitrary initial data A{0) G 
^-yiJ^u), it is a generalized solution. 
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The average values of observables (mean values of observables) are defined by the positive 
continuous linear functional on the space S,{J-'n) 

oo ^ 

{A){t) = {A{t),D{0)) = (/,D(0))-i5^-Tri,...,„A„(t)D°, (10) 

n=0 

where Tri_..._„ are the partial traces over 1, . . . , n particles, -D(O) = (1, D^, . . . , . . .) is a se- 
quence of self-adjoint positive density operators defined on the Fock space J^-^ that describes the 
states of a quantum system of a non-fixed number of particles and (/, -D(O)) = X^^^o^'^^i. - .^-^n 
is a normalizing factor (the grand canonical partition function). For -D(O) G £^(J-'-^) and 
A{t) G 'C(^-h) average value functional ( ITO!) exists and determines a duality between observ- 
ables and states. 

We note that in case of a system of fixed number of particles the observables and states are 
one-component sequences A^^\t) = (0, . . . , 0, AN{t), 0, . . .) and D^^\Q) = (0, . . . , 0, D%, 0, . . .), 
respectively, and therefore, average value formula ffTOj) reduces to the functional 

(A(^))(t) = {Ti,_^D%)-'Ti^_j,Aj,{t)D%. 

It is usually assumed that the normalizing condition: Tri^ jy-D^ = 1, take place. 

On the spaces ii('Hj) mapping ([7]) is defined and have the similar properties as in case of 
the Maxwell-Boltzmann statistics stated above. 

Symmetrization and anti-symmetrization operators ((21) are integrals of motion of the Heisen- 
berg equation (jl]), and as a consequence it holds 

and 

where the operator Afn is defined by (El). Thus, the symmetry of observables is preserved in 
process of the evolution. 

2.3 The von Neumann equation: the evolution of states 

As a consequence of the validity for functional (fTOl) of the following equahty 

oo ^ 

{A{t), D{0)) = (/, D(0))-i Tn,...,. Gnml Dl = 

n=0 

oo _ 

= (/, Gi-mo))-' Yl Tn,...,. < gn{-t)D'n ^ (A mr\A{o), D{t)), 

n=0 

where D{0) = (1, -D?(l), . . . , -D°(l, . . . ,n), . . .) G £,^{J^y), it is possible to describe the evolution 
within the framework of the evolution of states. Indeed, the evolution of all possible states, 
i.e. the sequence D{t) = (1, Di(t, 1), . . . , Dn{t, 1, . . . , tt,), . . .) G £,^{J^-u) of the density operators 
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Dn{t), n > 1, is described by the initial- value problem of a sequence of the von Neumann 
equations (the quantum Liouville equations) [3lll] 

j^D{t) = -^D{t), (11) 
D{t)\,=o = D{0). (12) 

The generator (—A/") = ©J^o(~-^n) ^^e von Neumann equation (fTTj) is the adjoint operator 
to generator of the Heisenberg equation (jl]) in the sense of functional f llUp . and for / G 
^liJ^n) C V{Af) C S^^i^n) it is defined by the formula 

(-Ar„)/„ = -t{Hr,fn - fnHn), (13) 

where the operator i/„ is the Hamiltonian ([3]). 

On the space of sequences of trace class operators £^(J^-^) for initial- value problem (fTTi) - (fT^ 
the following statement is true. 

Theorem 2. A unique solution of initial-value problem f|TTl) - f|T2l) of the von Neumann equation 
is determined by the formula 

D{t)=g{-t)D{0), (14) 

where the one-parameter family of operators g{—t) = Q)'^=QGn{—t) , is defined by 

6^„(-t)/n = e-**^"/„e**^". (15) 

If D{0) E S^KJ^'u) C £,^{J^-}{), operator (IT^ is a strong (classical) solution and for arbitrary 
-D(O) G £,^{J-''n) it is a weak (generalized) solution. 

On the space il^('H„) the mapping t — )■ Qn{—t)fn is an isometric strongly continuous group 
that preserves positivity and self-adjointness of operators. If /„ G £.l((Hn) C V{—Afn), in 
the sense of the norm convergence of the space il^(7{„) there exists the limit by which the 
infinitesimal generator of the group of evolution operators f[T^ is determined by 

lim i(^„(-t)/„ - /„) = -i{Hnfn - fnHn), 

where Hn is the Hamiltonian ([3]) and operator f[T^ is defined on the domain I){Hn) C Hn- 
On the spaces il^('H^) mapping (fT5|) is defined and have the similar properties as in case of 

the Maxwell-Boltzmann statistics stated above. 

Symmetrization and anti-symmetrization operators ([2]) are integrals of motion of the von 

Neumann equation (ITTl) . and as a consequence it holds 



g^{-t)s^ = s^gn{-t), 

and 

where the operator {—Mn) is defined by formula (fT3l) . Thus, the symmetry of states is preserved 
in process of the evolution. 
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3 The evolution equations for quantum correlations 

One more an equivalent approach to the description of the evolution of states of quantum many- 
particle systems is given by means of operators which are interpreted as correlation operators. 
The correlation operators are governed by the von Neumann [32], [41]. To justify such hierarchy 
of quantum evolution equations, at first we consider in detail the motivation of the description 
of states within the framework of correlation operators or in other words, we consider the origin 
of the microscopic description of correlations in quantum many-particle systems. 

3.1 Correlation operators 

We introduce a sequence of operators g{t) = (0, gi(t,l), . . . , gs{t, 1, . . . , s), . . .) G defined 
by the cluster expansions of the density operators D{t) = (1, 1), . . . , Dgit, 1, . . . , s), . . .) G 

D,{t,Y) = g,{t,Y) + J2 II 9\x,\it^X,), s > 1, (16) 

|P|>1 

where J2p-y=\j x, |P|>i possible partitions P of the set F = (1, . . . , s) into 

|P| > 1 nonempty mutually disjoint subsets Xi C Y. The operators defined by recursion 
relations f|T6l) are interpreted as the correlation operators of quantum many-particle systems, 
i.e. operators that characterize the correlations of particle states. 

In order to construct a solution of recursion relations flTB]) . i.e. to express the correlation 
operators injerms of the density operators, we introduce some notions. On sequences of 
operators /, / G £^(J-'-h) we define the *-product 

(/ * 7)|y|(r) = J2 MZ) f\Y\Zl{Y\Z), (17) 
ZCY 

where J2zcy subsets Z of the set 1" = (1, . . . , s). By means of definition (ITTl) 

of the *-product we introduce the mapping Exp^ and the inverse mapping Ln^, on sequences 
h = (0, ^i(l), . . . , /in(l, . . . , n), . . .) of operators /i„ G il^('H„) by the expansions 

(Exp, h)\y\{Y) = {l + J2—) = (18) 

71. 

n=l 

= /5|y|,0+ Yl n V.lTO. 
P-Y=[J,X, X,cP 

where we use the notations accepted in (fT6l) . S\y\,o is a Kronecker symbol, I = (1, 0, . . . , 0, . . .), 
and respectively, 

(Ln.(I + /.))|y|(F) = {J2^-ir~'—)\Yi(Y) = (19) 

ri=l 

= E (-i)'"-'(ipi-i)! n v.iM- 
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Hence in terms of sequences of operators recursion relations ( IT6l) are rewritten in the form 

D{t)=Exp,g{t), (20) 

where D(t) = I + {0, Di(t,l), . . . , Dn(t, 1, . . . ,n), . . .). From this equation we obtain 

g{t) = Ln, D{t). 

Thus, according to definition f|T71) of the *-product and mapping f|T9l) . in the component- wise 
form solutions of recursion relations (|T6l) are represented by the expansions 

gsit,Y) = Ds{t,Y) + Yl (-I)I^I-'(IPI-I)! n ^>1- (21) 

P:Y = \J^X„ X,CP 
|P|>1 

The structure of expansions fl2ip means that the correlation operators have a sense of cumulants 
(semi-invariants) of the density operators governed by the Cauchy problem of the von Neumann 
equations (fTT]) - (fT2|) . Therefore correlation operators (12T|) give an alternative approach to the 
description of the state evolution of quantum many-particle systems, namely within the frame- 
work of dynamics of correlations. 

We emphasize that the possibility of the description of states within the framework of corre- 
lations arises naturally as a result of dividing of the series in expression fllOj) by the normalizing 
factor series, i.e. in consequence of redefining of functional ffTOj) . 

To justify this proposition, i.e. to construct the mean- value functional in terms of correlation 
operators fl2Tl) . we introduce necessary notions and formulate some equalities. For arbitrary 
/ = (/o) /i) • • • ) /n, • • •) G il^(J-'-^) and y = (1, . . . , s) we define the linear mapping dy : f 
Dyf, by the formula 

(fy/)„ = fiY\+n{Y, |y| + 1, . . . , \Y\ +n), n> 0. (22) 
For the set {Y} consisting of one element Y = (1, . . . , s) we have respectively 

(f {y}/)n = fi+ni{Y}, s + l,...,s + n), n > 0. (23) 
On sequences Oy/ and dyf we introduce the *-product 

(Oy/*^yJ)|x|(X) = J2 f\Z\ + \Y\iY,Z)f^x\Z\ + \Y'\iY',X\Z), (24) 

ZCX 

where X, Y, Y' are the sets, which elements characterize clusters of particles, and YIzgx 

sum over all subsets Z of the set X. In particular case {Y = 0, y = 0) definition (!24|) reduces 

to ([nD. 

For / = (0, /i, . . . , /n, . . .), /„ G S^iHn), according to definitions of mappings (fTSjl and fl23|l . 
the following equality holds 

(){y}ExpJ = ExpJ*t)|y}/. (25) 
and for mapping fl22|) correspondingly, 

OyExpJ = ExpJ* fxj*...*c)x|p|/, (26) 

P: y=[J^ Xr 
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where J2p Y=\j-Xi ^^^^ possible partitions P of the set F = (1, . . . , s) into |P| 

nonempty mutually disjoint subsets Xi C Y. 
According to the definition 

oo ^ 

^-^ n\ 

n=0 

for sequences f,fE the equality holds 

il,f*f) = {lj){lj). (27) 

In terms of mappings (1221) and (1231) the generalized cluster expansions of solution (fT4|) of the 
von Neumann equation 

x\y)= 5^ n^(^'^*)' ^^1' (28) 

P:({y},x\y)=ux,x,cP 
where X \ 1" = (s + 1, . . . , s + n), take the form 

r)yZ^(t) = r){y}Exp, (7(t)- (29) 

We now deduce from the definition of functional flTUl) the functional for average values of 
observables in terms of correlation operators (I^Tj) . for example, for s-ary observables = 
(0, . . . , 0, a^(l, . . . , s), . . . , <...<i.=i '^^(^1' . . . , Zs), . . .), i-e. 

= i(/,D(t))-iTrya,(F)(/,r)yD(t)), 

where Try = Tri^... s. Using generalized cluster expansions ( 129|) and as a consequence of equal- 
ities (ESD and d^, we find 

{I,D{t))-\l,dYD{t)) = {I,D{t))-\l,d{y}Exp, git)) = 

= {I,D{t))~\l,ExpMt)*d{y}g{t)) = {I,D{t))-\l,Exp^g{t)){I,-0{y}g{t)). 

According to generalized cluster expansions fl29l) . as a final result we derive the following rep- 
resentation of the mean-value functional in case of s-ary observables 

(A(^))(t) = iTrya,(F)(/,c)|y|<7(t)), 

or in the componentwise form 

1 °° 1 

= - ^ - Tri„„,,+„ asiY)g,+4t, {¥}, s + 1, . . . , s + n), (30) 
s. n. 

n=0 

where the correlation operators of particle clusters gi+nit) are defined as solutions of generalized 
cluster expansions ( l29ll . namely 

g^Ut,{y},X\Y)= Yl (-I)'^I-^(IPI-I)! n ^(^'^^)' ^ > 1' (31) 

P:({y},x\y)=u,x. x,cp 

where the density operator D{t,Xi) is solution f lT4|) of the von Neumann equation (fTTI) . For 
A(^) e Il(JW) and e HH-^w) functional ([30D exists. 

It should be emphasized that correlation operators that belong to the spaces £^(J^-^) describe 
only finitely many particles, i.e. systems with finite average number of particles. 



Theory of quantum evolution equations 



12 



3.2 The von Neumann hierarchy 

We consider quantum systems of particles obeying the Maxwell-Boltzmann statistics with the 
Hamiltonian 

n n n 

H^ = Y,m + Y. E '^^'n^i,---,^.), (32) 

i=l k=l ii<...<ifc=l 

where ^'•'^^ is a A;-body interaction potential. The evolution of all possible states is described in 
terms of the correlation operators g{t) = (0, gi(t,l), . . . , gs{t, 1, . . . , s), . . .) G £}{J^y) governed 
by the following Cauchy problem 

^^g^{t,Y) = -K{Y)gs{t,Y) + (33) 

+ E E ••• E {-^f■F'^\z,,...,Z\^\))l[g\x.\{t,x,), 

P:Y = \J^Xi, ZiCXi, Z|p|CX|p|, XiCP 

|P| > 1 ^1 5^ Z|p| ^ 

where the operator (— A/'j^"'') is defined by 

J2 p ■ Y = \j-x^ \p\ > 1 possible partitions P of the set Y = (1, . . . , s) into 

|P| > 1 nonempty mutually disjoint subsets Xi C Y and J2z cx z 7^0 ^ over nonempty 
subsets C Xj. We refer to the hierarchy of equations (!33|) as the von Neumann hierarchy [32] . 

It should be noted that the von Neumann hierarchy (133|) is the evolution recurrence equations 
set. We cite an instance of the typical equations of hierarchy (l33l) 

j^g,{t,l) = -Af{l)g,it,l), 

j^g2it, 1, 2) = -Ar2(l, 2)(72(t, 1, 2) - Ar£)(l, 2)(7i(t, l)g,{t, 2), 

^^3(t, 1, 2, 3) = -A/'3(l, 2, 3)(73(t, 1, 2, 3) - 
at 

-{Af^^^{l,2)+Ar^^;{l,3)+Ar!^^il,2,3))g,{t,l)g,{t,2,3)~ 

-(Ar(?(l, 2) + Ar(f)(2, 3) + Ar(?(l, 2, 3))<7i(t, 2)g,{t, 1, 3) - 

3) + A/;i?(2, 3) + Ar£)(l, 2, 3))<7i(t, 3)(72(t, 1, 2) - 

-Ari2)(l, 2, 3)(1, 2, 3)(7i(t, l)(7i(t, 2)(7i(t, 3). 

In case of a two-body interaction potential ([3]) the von Neumann hierarchy fl33l) reduces to the 
following form 

^(/,(t,r) = -A/;(r)(7.(t,y)+ (36) 

+ E 5^ 5^ (-Ari,t(Zl,22))^7|Xi|(t,Xi)(7|X.|(t,X2), S>1, 
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where ^p.y=Xi u^2 possible partitions P of the set y = (1, . . . , s) into two 

nonempty mutually disjoint subsets Xi G Y and X2 C Y. In terms of kernels of correlation 
operators the first two equations of the von Neumann hierarchy (136|) have the form 

d 1 

d . , 
i-Q^92{t,quq2;qi,q2) = 

1 ^ 

- 92) - - q2))gi{t, qu q'i)9i{t, gs; ga)- 

We remark that the von Neumann hierarchy (l33l) (or (|36l) ) is rigorously derived on basis of 
solutions fl2Tl) of cluster expansions f|T6|) of density operators governed by the von Neumann 
equation ffTTl) . 

To construct a solution of the Cauchy problem fl33|) - fl34|) we first consider its structure for 
physically motivated example of initial data, namely initial data satisfying a chaos property. A 
chaos property means the absence of state correlations in a system at the initial time. In this 
case the sequence of initial correlation operators is the following one-component sequence 

(7(0) = (0,(7?(1),0,...). (37) 

In fact, in terms of the sequence of the density operators this condition means that -D(O) = 
(l,D5'(l),L'5'(l)L'?(2),...,n"^^D°(i),...) in case of the Maxwell-Boltzmann statistics, then 
from ( j2TI) we derive ( !37ll . 

On basis of representation (I2T1) of the correlation operators in terms of the density operators 
and formula ([T3j) we have 

gs{t,Y)= (-i)"''"'(ipi-i)! n ^i^-i(-^'^^)ri^?(^)' 

P:Y = \J^X, X,CP i=l 

Taking into account the equality: -D?(^) = fi'i(^), 1 < i < s, we derive the formula of a solution 
of the Cauchy problem (1331) - (|34l) for initial data (1371) 

g,{t,Y)=Ql,{-t,Y)l[g',{t), s>l, (38) 

i=l 

where 2ls(— t, Y) is the sth-order cumulant of the groups of operators (ITSll of the von Neumann 
equations (ITTll 

2t.(-t,y)= Yl (-I)l'l-^(IPI-I)! n ^l^^lH'^^)' (39) 

P:Y=\J,X, X,CP 

and we use accepted above notations. It should be emphasized that the structure of expan- 
sions (!39l) of the evolution operators 2ts(t,y) means that they are determined by the cluster 
expansions of the groups of operators (fT5|l of the von Neumann equation (fTTjl . 
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Thus, the cumulant nature of correlation operators induces the cumulant structure of a one- 
parametric mapping generated by solution ( 138|) . From ( 138|) it is clear that in case of absence 
of correlations in a system at initial instant the correlations generated by the dynamics of a 
system are completely governed by cumulants (!39|) of groups (fTSjl . 

Let us indicate some properties of a cumulant (semi-invariant) (139|) of groups of operators 
(IT3]) . The corresponding-order cumulant is a solution of the cluster expansions of groups of 
operators Qs{—t) (HM 

gsi-t,Y)= Yl n ^^1' (40) 

P:y=U,X, X,CP 

where is the sum over all possible partitions P of the set F = (1, . . . , s) into |P| nonempty 
mutually disjoint subsets Xj C Y. The simplest examples of equations from recurrence relations 
fHOj) are given by the expressions 

6;i(-t,i) = 2ti(-t,i), 

g2i~t, 1, 2) = 2l2(-t, 1, 2) + 2li(-t, l)2li(-t, 2). 

In case of a quantum system of non-interacting particles we have: 2ts(— t) = 0, s > 2. Indeed, 
for non-interacting quantum particles it holds 

s 
i=l 

and hence, 

P:Y=\J.Xi XiCPli=l 
s s 

= J2{-i)'-Hs, k){k-i)\l[ g,{-t, i) = 0. 

fc=l i=l 

Here s(s, k) are the Stirling numbers of the second kind and the following equality is used 

s 

E (-l)"''"'(|P| - = ^)(^ - = '^M, (41) 

where 5^,1 is a Kroneker symbol. 

If s = 1, for /i G -Co(^) -C^('H) in the sense of the norm convergence of the space 
the generator of first-order cumulant (15^ is given by operator ([T5|l 

lirnl(2li(-t,l)-/)/i = -M/i. 
For s > 2 we obtain the following equality in the same sense 

hni \ 2l,(-t, 1, . . . , s)/, = (1, . . . , 5)/„ 
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where the operator (— A/'^t'') is defined by formula (!35|) . 

If fs G S^iTis), then for sth-order cumulant (!39|) of groups of operators (fT5|) the estimate is 
valid 

l|2i.H)/.L.(«.) < E m-m\fsLHn.)< (42) 



fc=i 



where s(s, k) are the Stirling numbers of the second kind. 

Hereafter we use the following notations: ({Xi}, . . . , {X|p|}) is a set, elements of which 

are |P| mutually disjoint subsets C F = (1, . . . , s) of the partition F : Y = U^^^Xi, i.e. 
|({Xi}, . . . , {X|p|})| = |P|. In view of these notations we state that {Y} is the set consisting 
of one element Y = of the partition P (|P| = 1) and \{Y}\ = 1. We define the 

declasterization mapping 9 : ({^i}, . . . , {X|p|}) — Y, by the formula 

^({Xi},...,{X|p|}) = y. (43) 

For arbitrary initial data a solution of the Cauchy problem ( p3|) -f p4|) of the von Neumann 
hierarchy is given by the expansion 



gs{t,Y)= 2l|p|M,{^i},---,{^|P|}) n ^^1- (44) 

P:Y=\J^X, X,CP 

Here 2l|p|(— t) is the |P|t/i-order cumulant of groups of operators f|T5|) defined by the formula 

2l|p|(-t,{Xi},...,{X|p|})= (45) 
E (-I)I^'I-^(IP'I-I)! n G\eiz,)ii-t,9iZ,)), 

P':({Xi},...,{X|p|})=Ufc^fc ZkCP' 

where ^p'.({Xi} {X|p|})=UfcZft possible partitions P' of the set ({Xi}, . . . , 

{X|p|}) into |P'| nonempty mutually disjoint subsets C ({Xi}, . . . , {X|p|}). The simplest 
examples of correlation operators (jUj) are given by the expressions 

^i(t,l)=2ti(-t,l)(7?(l), 

g2{t, 1, 2) = 2ti(-t, {1, 2})(7°(1, 2) + 2l2(-t, 1, 2)(7?(1)^?(2), 
g^it, 1, 2, 3) = 2li(-t, {1, 2, 3})^3°(1, 2, 3) + 2l2(-t, {2, 3}, 1)^?(1)(72°(2, 3) + 
{1,3}, 2)g%2)g',{l, 3) + 2t2(-t, {1,2}, 3)^?(3)^2°(1, 2) + 
+2l3(-t,l,2,3)^7?(l)^7?(2)^7?(3). 

The validity of solution expansion f l44|) can be verified by straightforward differentiation by 
time variable and also in the following way. Taking into account the fact that the von Neumann 
hierarchy fl55]) is the evolution recurrence equations set, we can construct a solution of initial- 
value problem ( I33p -f l5^ by integrating each equation of the hierarchy as the inhomogeneous 
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von Neumann equation. For example, using formula ( ITSl) . as a result of the integration of the 
first two equations of hierarchy ( l33l) we obtain the following equalities 

(7i(t,l) = ^i(-t,l)^°(l), 

^2(t,l,2) = 6;2(-t,l,2)(?°(l,2) + 

t 

+ j dhG2{-t + hA,2){-UU^,2))Gi{-h,l)Gi{-h,2)gl{l)gl{2). 



Then for the second term on the right-hand side of this equation an analog of the Duhamel 
equation holds 

t 

j dtiG2{-t + ti,l,2){-MU^,2))gi{-ti,l)Gi{-h,2)= (46) 



t 

/I 
dti—{g2{ti, 1, 2)g,{-h, i)g,{-ti,2)) = 



= g^i-t, 1,2)- g^{-t, i)gi{-t, 2) = ^^{-t, 1, 2), 

where 2l2(— t) is the second-order cumulant of groups of operators ffT^ defined by formula fl5I?]) . 
For s > 2 a solution of the Cauchy problem (l33l) - (!MI) constructed by iterations is represented 
by expansions (H4|) as a consequence of transformations similar to an analog of the Duhamel 
equation (l46ll . 

We note, that in case of initial data ( |371) solution (l38l) of the Cauchy problem (l33l) - (!Ml) of 
the von Neumann hierarchy may be rewritten in another representation. For tt, = 1, we have 

g,{t,l) = Ql,{-t,l)g'({l). 

Then, within the context of the definition of the first-order cumulant, t), and the dual 
group of operators ^i{t), we express the correlation operators gs(t), s > 2, in terms of the 
one-particle correlation operator gi{t) using formula f l38p . Hence for s > 2 formula (138 p is 
represented in the form of the functional with respect to one-particle correlation operators 

s 

gs{t, Y I gi{t)) = 2t,(t, Y) J] g^{t, i), s > 2, 

i=l 

where 21^(^,1") is st/i-order cumulant (1390 of the scattering operators 

s 

g^it, Y) = g^i-t, Y) II g^it, z), s > i. (47) 

i=l 

The generator of the scattering operator gt{Y) is determined by the operator 

, s s 

k=2 ji<...<ifc=l 
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where the operator (— A/'int^) acts on S^KHs) C Sl^^Hs) according to formula (135|) . 

On the space for initial- value problem ( 133|) - (!34|) of the von Neumann hierarchy the 

following statement is true [32], |41] . 

Theorem 3. For t G M a solution of the Cauchy problem fl33p -flM |) is given by the following 
expansion 

g,{t,Y) = g{t;Y\g{0))= (48) 
= {^i}' • • • ' n 9^x.\m, s > 1, 

where 2l|p|(— is the \P\th-order cumulant f l45|) of the groups of operators. For G 2j\{'Hn) C 
n > 1, expansion fHS]) zs a strong (classical) solution and for arbitrary initial data 
g^ G ii^ n> 1, it is a weak (generalized) solution. 

The intrinsic properties of constructed solution (HHj) are generated by the properties of cumu- 
lants (H5|) of groups of operators of the von Neumann equations. We indicate some properties 
of the one-parameter mapping: t — )■ Q{t\f), generated by solution psj) . 

For fs G £^('Hs), s > 1, the mapping Q(t; Y\f) is defined and, according to the inequality 



||2l|p|(-t, {X^}, {X|p|})/.||^,(^^^ < |P|! el^i||,.||,,(^^), 
the following estimate is true 

\\^it;Y\f)\\^.^^^^<s\e'^c% (49) 

where c = e^max(l, maxp.y=u^Xi ||£i(«|;^,|))- The mapping Q{t\f) has the group property, 
i.e. for / G £^(J-'-^) the equalities are fulfilled 

Gih + t2\f) = g{ti\g{t2\f)) = g{t2\g{ti\f)). 

On the subspaces ^liTis) C S^iTig), s > 1, the infinitesimal generator Af(Y\f) of the group 
Q{t; Y\f) is defined by the operator 

Ar(F|/) = -Ar,(y)/,(F) + (50) 

p 

+ E E ••• E {-MF'^\zi,---,z\^\))i[f\x^\{x,), 

P:Y = [J^X,, ZiCXu Z|p|CX|p|, X,CP 

|P| 7^ 1 Zi ^ Z|p| ^ 

We remark that a particular solution of the steady von Neumann hierarchy is a sequence of 
the Ursell operators, for example, its first two elements have the form 



9iit, 1) = e 



-I3K{1) 



^2(M,2) = e"'£''^^e-/^*(^-^)-/), 
where (3 is a parameter inversely proportional to temperature. 
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For the purpose of further apphcation we introduce the notion of correlation operators of 
clusters of particles 

g^+^{t, {Y},X\Y)= g{t- {Y},X\ F |(^(0)) = (51) 

= E 2i|Pi(-t,mxo},...,mxiP|)}) n*.iTO, 

P:({y},x\y)=uxi x,cp 

where 2l|p|(— t) is the |P|t/i-order cumulant defined by formula (HSj) . 

We remark that the relations between correlation operators of particle clusters '^)^^^^g{t) G 
£^(©J^o^s+n) correlation operators of particles fHHl) are given by the following equahties 

3i+n(t,{y],X\Y)= (52) 

= (-i)"-i-(iPi-i)! n E n si-,.i(*>2..)- 

P:({y},x\y)=ux, X,CP P':e(X,)=U^^ z,^ Z,,CP' 

In particular case n = 0, i.e. the correlation operator of a cluster of |y| particles, these relations 
take the form 

P:y=U,X,X,CP 

Due to cluster expansions fll6p it is possible to express many-particle correlation operators 
through the two-particle and one-particle correlation operators 

^i+„(t,{r},x\F) = 

p 

ZcX\Y,Zj^tD P:X\{YUZ)=\J.Xi i=l 

We note also that in case of many-particle systems obeying quantum statistics, i.e. many- 
particle systems of fermions or bosons, the von Neumann hierarchy has the form 

j^g,{t,Y)=Ar{Y\g{t)), s > 1, (53) 

where the hierarchy generator is determined by 
AfiY\g) = -Ms{Y)gs{Y) + 

+ E E ••• E {-^F'^\zu---.Z\^\))Sf\{g\xA{X;), 

P:y = U,^>, ^iCXi, Z|p|CX|p|, X,CP 

IPIt^I ^17^0 2|P|7^0 

and the operators Sf are defined by formula ([2]). A solution of the Cauchy problem of hierarchy 
(!53|) is given by the corresponding expansion 

gs{t,Y)= E 2l|p|(-t,{Xi},...,{X|p|})5f n ^^Ix.lM, ^>1- 

V■.Y={j^X, X,cP 
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4 The hierarchies of quantum evolution equations 

For the description of the evolution of infinite-particle quantum systems the hierarchies of evo- 
lution equations for marginal observables, marginal density operators and marginal correlation 
operators are used P, [3]. They are constructed as the evolution equations for one more but 
an equivalent method of the description of states and observables of finitely many particles. 

Usually the evolution of infinite-particle quantum systems is described within the framework 
of the evolution of states by the quantum BBGKY hierarchy for marginal density operators. 
An alternative method of the description of state evolution is given in terms of the nonlinear 
quantum BBGKY hierarchy for marginal correlation operators. An equivalent approach to the 
description of the evolution of quantum systems is given within the framework of the marginal 
observables governed by the dual quantum BBGKY hierarchy [21], [40j . 

For a system of a finite average number of particles there exists an equivalent possibility to 
describe observables and states, namely, by sequences of marginal observables (the so-called s- 
particle observables) B(t) = {Bq, Si (t, 1), ... , Bs{t, 1, . . . , s), . . .) and marginal states (s-particle 
density operators) -F(O) = (1, F^{1), . . . , -F°(l, . . . ,s), . . .) [T], [1]. These sequences are corre- 
spondingly introduced instead of sequences of observables A{t) and density operators -D(O), in 
such way that mean value (fTOjl does not change, i.e. 

CO ^ 

{A){t) = {I,Dm-\A{t),Dm = (/,D(0))-i^-Tri„„,„A„(t)D° = (54) 

n=0 

oo _ 

= {B{t),F{0)) = 5^ - Tri,„„, B,{t, 1, . . . , 5) F°(l, . . . , 5), 

where (/, D(0)) = ^^o^'^^iv.n-^n a normalizing factor and / is the identity operator. 
Thus, the relationship of marginal observables and observables is determined by the formula 

B,{t,Y) = y^^-^ V A.n{t,Y\Ui,...,jn)), s>l, (55) 
^ — ^ n! ^ — ^ 

n=0 ji^...^j„ = l 

where Y = {1, . . . ,s). In terms of the density operators the marginal density operators are 
defined respectively 

00 ^ 

(y) = (/, D(0))-i J2 -J^i,-,s+n I^°+.(l, ...,s + n), s>l. (56) 

n=0 

As we can see from (l55l) one component sequences of marginal observables correspond to observ- 
ables of certain structure, namely the marginal observable B^^^ = (0,ai(l),0, . . .) corresponds 
to the additive-type observable A^^'> = (0, ai(l), . . . , ^^=1^1 (^), • • ■), and in the general case the 
/c-ary-type marginal observable B^''^ = (0, . . . , 0, afc(l, . . . , fc), 0, . . .) corresponds to the fc-ary- 
type observable A^''^ = (0, . . . , 0, 0^(1, . . . , /c), . . . , Y.Z<-<ik=i ■ ■ ■ ^^k), ■ ■ ■)■ 

We emphasize that the evolution of marginal observables flS5]) of both finitely and infinitely 
many quantum particles is described by the initial- value problem of the dual BBGKY hierarchy. 
For finitely many particles the dual quantum BBGKY hierarchy is equivalent to the Heisenberg 
equations (j4]) (the dual equation to the Heisenberg equation (jl]) is the von Neumann equation 
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4.1 The dual quantum BBGKY hierarchy for marginal observables 

The evolution of marginal observables is described by the initial-value problem of the dual 
quantum BBGKY hierarchy 

, s s 

-B,{t,Y) = + E Mnt{ji,j2))B.,{t,Y)+ (57) 

i=i ii<i2=i 

s 

+ Yl -Mnt(jl,j2)i?s-l(t,n(jl)), 
jl¥=j2 = l 

Bsit) |i=o= Bl s > 1, (58) 

where on ^oiTin) C 2{'Hn) the operators and Mmtiji,j2) are correspondingly defined by 

formulas 

W)9n = -i{gnK{j) - K{j)gn), (59) 

We refer to recurrence evolution equations (1571) as the dual quantum BBGKY hierarchy since 
it is the adjoint hierarchy of evolution equations to the quantum BBGKY hierarchy for the 
marginal density operators [3., 4J with respect to bilinear form f fTOj) . In case of the space H = 
L'^{W^), evolution equations fl57|l for kernels of the operators Bs{t), > 1, are given by 

t ^^B,{t, gi; q[) = -^(-A,, + A,. )Si(t, q,; q[), 



. d 



1 

Bs{t, gi, g.; gi, ...,?;) = (- - $^(-A,^. + A,. ) + 



+ J2 i^i<l'h-<lh)-'^(lh-<ln)))Bsit,qi,...,qs;q[,...,qi) + 

i=ii<i2 

s 

i=ii7^i2 

We note that in case of many-body interaction potentials f l32p hierarchy fl57p has the form 

^Bs{t,Y)=K{Y)Bs{t,Y) + 
at 

+ E 713^ E ^f^nhjl,---,J^)Bs-n{t,Y\{n,...,J^}), S>1. 

n=l k=n+l ^ ' j^^...^jk=l 

On space 2,(l-ik) the operator M^^^ is defined by 

M^'hk = -^{9k^^'^-^^'^9k). (60) 
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We note that for finitely many particles the dual quantum BBGKY hierarchy ( 1571) (or ( l60l) ) is 
rigorously derived on basis of the Heisenberg equation (jll) according to definition (1551) . Hence 
the structure of evolution equations for marginal observables is defined by the structure of 
expansion (!55|) . 

Let us introduce some abridged notations: Y = {1, . . . , s), X = {ji, . . . , j„) C Y and {Y\X} 
is the set consisting of one element Y\X = (1, . . . , s)\(ji, . . . , i.e. the set {Y\X} is a 
connected subset of the set Y. 

To construct a solution of the Cauchy problem (I57|l -f l58|) we introduce the (1 + n)t/i-order 
cumulant of groups of operators ([7]) as follows [30] , [33] 

^,+n{t,{Y\X},X)= (61) 
E (-I)I^I-^(IPI-I)! n ^^0' 

P:{{Y\X},X)=[J^X, X,CP 

where is the sum over all possible partitions P of the set {{Y\X},ji, . . . ,jn) into |P| 
nonempty mutually disjoint subsets Xi C {{Y\X},X). For example, 

^,{t,{Y}) = gs{t,Y), 

Mt,{y\U)},j) = Os{t,Y) - g,_,{t,Y\{j))g,{t,j). 

Let us indicate some properties of cumulants fl^ . If n = 0, for G £o(^s) C il('Hs) in the 
sense of the *-weak convergence of the space 2{T-Ls) the generator of first-order cumulant ( ETj) 
is given by operator ([6]) 

lim ^(2li(t, {Y}) - I)gs{Y) = KgsiY). 

In case of n = 1 for gg G £o('H<i) C 2>{1-Ls) we obtain the following equality in the sense of the 
*-weak convergence of the space S.{l-Ls) 

\im-^Mt,{Y\{j)},j)g,{Y)= J] UUh3)9s{Y), 

ie(n(j)) 

where the operator Mint{i,j) is defined by f l59l) . and for n > 2 as a consequence of the fact that 
we consider a system of particles interacting by a two-body potential ([3]), it holds 

lim 1 2ti+„(t, {Y\X}, X)g,{Y) = 0. 



Correspondingly in case of n-body interaction potential fl32|) it holds 

Jim 7 2l„(t, 1, . . . , n)gn = J^tt 9n, 

where the operator M!^^^ is defined by formula (!60l) . 

If gs € £(7^<j), then for (1 -|- n)t/i-order cumulant (16T]) of groups of operators ([7]) the estimate 
is valid 

||2li+„(%.||^(^^) < (|P|-l)'lkL(H.)^ (62) 

P:({y\X},X)=UX. 

71+1 

< 5^ s(n + 1, A:)(A; - 1)!||^7.||,(^^) < ^!e"+'|k.„,(^^), 



fc=i 
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where s{n + l,k) are the Stirhng numbers of the second kind. 

On the space for abstract initial- value problem (l571)-(l58l) the following statement is 

true [10]. 

Theorem 4. 7/5(0) G and 7 < , then for t a unique solution of the Cauchy 

problem fl57|) - fl58l) of the dual quantum BBGKY hierarchy exists and it is determined by the 
expansion 

S -. s 

Bs{t,Y) = Y^- J2 ^i+nitAY\X},X)Bl^iY\X), s>l, (63) 

where the (1 + n)th-order cumulant 2li+„(t, {Y\X}, X) is defined by formula (1611) . For -B(O) G 
£°(J-'^) C £.j{J-'y) it is a classical solution and for arbitrary initial data B{0) G £^(J-'-^) it is a 
generalized solution. 

The simplest examples of marginal observables ( 163|1 are given by the expressions 
Si(t,l)=2li(t,l)i??(l), 

B2{t, 1, 2) = 2ti(t, {1, 2})i?0(l, 2) + 2l2(t, 1, 2)(i?0(l) + S?(2)). 

We remark that at the initial time t = solution fl63|) satisfies initial condition f l58|) . Indeed, 
according to definition ([7]) and equality fHT]) for tt, > 0, we have 

2ti+„(0, {n^},X) = 5^ (-I)I^I-'(IPI - 1)!/ = IKo- 

P:{{Y\X},X)=U,X, 

On the space £^(J-V^) a solution of the initial- value problem of the dual BBGKY hierarchy 
defines a one-parametric mapping with the following properties j33]. 

Owing to estimate ( 162|) for cumulants ( 16T]) . under the condition that 7 < we have 

II^WL.(^.)<^'(l-7e)-iS(0)||,^(^^). (64) 

Then on the space g G £,^{J^y) the one-parametric mapping 

R^3t^ {U{t)g)s{Y) = (65) 

= 5^7 TT V 2ti+„(t,{n^},x)^?,_„(y\^), s>i, 

n=0 ^ ^ ji^...^j,_„=l 

is a Cg-group. The infinitesimal generator P = ©^q-^s '^^ ^^^^ group of operators is a closed 
operator for the *-weak topology and on the domain of the definition T>{'D) C 2,^{J^T-i) which 
is the everywhere dense set for the *-weak topology of the space £^(J-'-^) it is defined by the 
operator 

{Vg),{Y)=Us{Y)g,{Y) + 

S s ^ s 

-^H-lH (k-n)\ ^ M^^i{jl,...,]k)gs-niY\{ji,...,Jn}), 
n=\ ■ k=n+l ^ ji^...^jfc = l 
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where the operator Af-j'^^ is given by formula (160|) . 

In capacity of initial data we consider the additive-type observables, i.e. the one-component 
sequences B^^^O) = (0,ai(l),0, . . .) (the k-aij marginal observable is the sequence B^''\0) = 
(0, . . . , 0, afc(l, . . . , A;), 0, . . .) [22]). In this case solution expansion attains the form 

s 

Bi'\t, y) = 2l,(t, 1, . . . , 5) J] ai(j), s > 1. (66) 

i=i 

For such additive-type observable as number of particles, i.e. one-component sequence A^(0) = 
(0, /, 0, . . .), according to definition of cumulants (EI]), solution expansion f l^ get the following 
form 

s 

{N{t))s{Y) = 2t,(t, 1, . . . , s) 5^ / = s 6s,il, s>l, 

where / is the identity operator and 6s,i is a Kroneker symbol. Hence we have 

|(iV(t),F(0))| = |TriF{'(l)| < llFfll^,^^^ < oo. 

Thus, the marginal density operators from the space £^(J-'-^) describe quantum systems of 
finitely many particles. 

We note that solution expansion fl63l) . i.e. nonperturbative solution expansion, can be derived 
from solution (Q of the initial- value problem of the Heisenberg equation (ll])-([5]) on the basis of 
expansions (|55|) . Since hierarchy fl571) has the structure of recurrence equations, we also deduce 
that the solution can be constructed by successive integration of the inhomogeneous Heisenberg 
equations. Solution (jHSD is represented in the form of the perturbation (iteration) series as a 
result of applying of analogs of the Duhamel equation to cumulants (I6T]) of groups of operators 

m- 

Cluster expansions (HOj) of group of operators can be put at the basis of all possible solu- 
tion representations of the dual quantum BBGKY hierarchy (1571) . In fact, solving recurrence 
relations (1401) with respect to the Ist-order cumulants for the separation terms, which are 
independent from the variable Y\X = {ji, . . . ,js-n) 



Ql,UtAY\X},X)=J2Mt,{Y\XUZ}) Yl (-1)"''|P|! 

ZCX P:X\Z=[J^X, i=l 



2li(t,{X,}), 

where Ylzcx ^ sum over all subsets Z G X of the set X and taking into account the identity 



J2 {-lf^\P\\ll'^iit,{X,})g.^^{Y\X)= Yl {-lf^n9s-n{Y\X), (67) 

P:X\Z=[J^X, i=l P:X\Z=\J^X, 



and the equality 



Y (-1)'" |P|! = (-1)'^^^', (6^ 

P:X\Z=\J^X, 
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for expansion ( l63i) we derive 

Introducing the operator o"*" (an analog of the creation operator [40j): 

{^^9)s{Y) = 9s-i{Y\{j)) (70) 

defined on 2^{J^y), as a result of the symmetry property of the Maxwell-Boltzmann statistics, 
expression (169|) can be rewritten in the following compact form 

OO ^ 71 I 

^(^) = E E (-1)""' k\(n-k)\ i''^T~'^(t){a^)'B' = e-^"Q{tV"B'- 

n=0 ' k=0 

We can obtain one more representation for a solution of the initial- value problem of the dual 
quantum BBGKY hierarchy, if we express the cumulants 2li+„(f), n> 1, of groups of operators 
([7]) with respect to the Ist-order and 2n(i-order cumulants. In fact, it holds 

|P| 

2ii+.(t,{(n^)},x)= Y.^,{t,{Y\x},{z}) {-if^nW'^^it.ix^M 

ZCX, V:X\Z=[j.Xi i=l 

where ^zcx, is a sum over all nonempty subsets Z C X of the set X. Then taking into account 

identity (1^ and equality (|^ . we get the following representation for solution expansion (1^ 
of the dual quantum BBGKY hierarchy 

5,(t,F)=2li(t,y)5°(F) + 

+E T^z-^x E E (-1)""^^' 2t2(t,{n^},{^})5L„(n^), 

Zj^d) 

where y = (1, . . . , s), X = Y\{ji, . . .Js-n), i-e. Y\X = {ji, . . .,js-n)- 
4.2 Marginal density operators 

As stated above f l54|) the mean value of the marginal observable B[t) e Sl^{J^y) at t G M in 
the initial marginal state F(0) = (1, F°(l), . . . , F°(l, . . . , s), . . .) e Ili(JW) is defined by the 
functional 

oo ^ 

(i?(t), F(0)) = 5^ - Tn,...,, 5,(t, 1, . . . , s)F^{l, ...,s). (71) 
.=0 ^■ 

According to estimate (El]), functional fl7T]) exists under the condition that 7 < e~^, and the 
following estimate holds 

|(S(t),F(0))|<e2(l-7e)-iS(0)||,^(^^)||F(0)||,,_^(^^^. 
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In consequence of the vahdity for functional ( ITTi) of the following equalities 

{I,Dm-\A{t),Dm = {I,D{t))-\A{0),D{t)) = (S(0),F(t)) = 

oo ^ 

= 5^ -Tri,...,sS°(l,...,5)F,(t,l,...,5), 

where the marginal density operator Fs{t,l, . . . , s) is defined by means of density operators 
(fT4|) according to formula (156!) and the marginal observable Bg{l, . . . ,s) is defined by means 
of ([5]) according to formula (155!) respectively, it is possible to describe the evolution within the 
framework of the marginal state evolution. 

According to fl30p . marginal density operators (1561) can be defined in terms of the correlation 
operators of clusters of particles by the expansion 

oo ^ 

{t,{Y},s + l,...,s + n), s>l, (72) 

n=0 

where the correlation operator 5'i+„(t) is given by expansion (l3Ti) . Starting from an alternative 
approach to the description of states by the von Neumann hierarchy, we define the marginal 
density operators by the use of solutions of the Cauchy problem of the von Neumann hierarchy 
for correlation operators of clusters of particles (151 p . As is obvious from what will be said 
the cumulant structure of the von Neumann hierarchy solution flSTl) induces the cumulant 
structure of solution expansion f l97|) of initial- value problem of the quantum BBGKY hierarchy 
for marginal density operators that are adopted for the description of infinite-particle systems, 
i.e. infinite-particle dynamics is generated by the dynamics of correlations. 

4.3 The quantum BBGKY hierarchy 

The evolution of states is described by the sequences F{t) = (1, Fi(t, 1), . . . , Fs{t, 1, . . . , s), . . .) 
of the marginal density operators that satisfy the Cauchy problem of the quantum BBGKY 
hierarchy 

j/s{t,Y) = (Y^i-mj)) + E (--MntO-i,j2)))F.(t,i^) + 

i=i ji<j2=i 

s 

+ J2 Trs+i{-Mnt{3, s + l))F,+i(t, Y,s + 1), (73) 

F,{t,Y)\t=o=F^{Y), s>l, (74) 

where Y = (1, . . . , s) and on QK'Hn) C Q^iTin) the operators {—Af{j)) and (— A/'int(ji, ^2)) are 
correspondingly defined by formulas 



(-A/'(j))/n = -^(^(j)/«-/n^(j)), 
{-Mnt){jl,32)fn = -«($(jl, j2)/n " fn^{jl,j2)). 



(75) 
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In case of the space "H = L'^{W^), hierarchy of evolution equations (173|) for kernels of the 
operators Fs{t), s > 1, (the marginal or s-particle density matrix) are given by 

1=1 

s 

+ {^((li-^J)-^((l'^-(l'j)))w,ql,■■■.qs;q'l,■■■,qs) + 

i<j=i 

+ Y dqs+i{^{qi - qs+i) - ^{q- - q^+i)) Fs+i{t, qi, . . . , qs, g.+i; q[,...,q'„ qs+i). 

We note that in case of many-body interaction potentials ( 132|) hierarchy ( 1731) has the form [33] 
j^F,{t, Y) = -^f,{Y)F,{t, Y) + (76) 
+ E;^Tr,+i,„.,,+„ 5^ (-Ar(lfl+"))(Z,s+l,...,s + n)F,+„(t), s>l, 

n=l ■ Z CY, 

where the operator M!^^^ is defined on S^KHn) C £,^{'Hn) by formula ( |60ll . 

For finitely many particles the quantum BBGKY hierarchy (|73l) is rigorously derived on basis 
of the von Neumann equation f|TT]) according to the definition of marginal density operators fl56|) . 
The rigorous derivation of the quantum BBGKY hierarchy from the von Neumann hierarchy 
is given in [16] . Another way of looking to the derivation of the quantum BBGKY hierarchy 
consists in the construction of the adjoint (dual) equations to the dual quantum BBGKY 
hierarchy flFTI) with respect to bilinear form flTTl) . 

On the space £^(J-'-^) for abstract initial- value problem (173l) - (17ill the following statement 
holds [31] • 

Theorem 5. If F{0) G -^^(J-'^) and a > e, then for t G M a unique solution of the Cauchy 
problem (!73l) - (!7il) of the quantum BBGKY hierarchy exists and is given by the expansion 

CO 

F,{t,Y) = ^-Tr,+i„„,,+„2li+„(-t,{y}, X\Y)Fl^{X), s > 1, (77) 

n=0 

where Y = {1, . . . ,s), X = (1, . . . , s + n), and the evolution operator 

Ql,U-t,{Y}, X\Y) = Y (-1)'""'(|P| - 1)! n O\eixd-t,0{X,)) (78) 

P :{{¥}, X\Y)=[j,X, X.CP 

is the (1 + n)th-order cumulant of groups of operators f[T^ . is the sum over all possible 
partitions P of the set {{Y}, X\Y) into |P| nonempty mutually disjoint subsets Xi C ({^}, X\ 
Y). For initial data F(0) G H^q ^ ^ii^u) it 

is a strong solution and for arbitrary initial data 

from the space £^(J-'-h) it is a weak solution. 
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Owing to estimate ( H2l) . i.e. 

for F° G series fITTI) converges on the norm of the space provided that a > e, 

and the inequahty holds 

||F(t)||£i(^«)<cJ|F(0)||.i(^„), 

where Cq, = e^(l — -)^^. The parameter a can be interpreted as the value inverse to the average 
number of particles. 

Let us indicate some properties of cumulants (|75i) . If n = 0, for fs G ^liTis) C ii^('Hs) in 
the sense of the norm convergence of the space -C^('Hs) the generator of first-order cumulant 
fl78|) is given by operator f|T3|) 

limi(2ti(-t, {y}) - I)UY) = -KfsiY). 

t^o t 

In case of n = 1 for Z^+i G £,l{'Hs+i) C £^(7{s+i) we obtain the following equality in the same 
sense 

hni 1 2t2(-t, {Y},s+ s + 1) = V(-A/'int(^, s + s + 1), (79) 

t^o t ^-^ 

i=l 

and for n > 2 as a consequence of the fact that we consider a system of particles interacting by 
a two-body potential ([3]), it holds 

hm - {F}, = 0. 

t->o t 

With a view to generality we consider the case of many-body interaction potentials fl32|) . On 
the space Sj]^{J^u) ^ solution of the initial- value problem of the quantum BBGKY hierarchy is 
defined a one-parametric mapping (the adjoint mapping to mapping flB5]) in the sense of bilinear 
form (HD) 

R'3t^ {U{-t)f),{Y) = V-Tr,+i,....,+„2li+„(-t,{F},X\y)/,+„,(X) 

n=0 

with the following properties [33]. If / G £^(J-'-^) and a > e, then one-parametric mapping (IHOll 
is a Co-group [33]. On the subspace S^^o C. £,l^{J^-u) the infinitesimal generator B = ^'^^^Bn 
of this group is defined by the operator 

{Bf),{Y) = -^,{Y)UY) + 

+ E E^Tr,+i,...,,+„(-Ar4f"))(zi,...,z,,X\y)/,^^ .>!, 

fc=l ■ ii^...^ifc=l n=l 

where on QKHs+n) C £^('Hs+n) the operator (— A/'int^"'') is defined by formula (l35!l . 
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We indicate that nonperturbative solution ( 177|) of the quantum BBGKY hierarchy is trans- 
formed to the form of the perturbation (iteration) series as a result of applying of analogs of the 
Duhamel equation ^ to cumulants (ITHj) of groups of operators. To reduce expansion (1771) to 
the BBGKY hierarchy iteration series we put groups of operators in the expression of cumulant 
(!78|) into a new order with respect to the groups of operators which act on the variables Y 

2ti+„(-t,{F},X\y)= (80) 

= J] %uzi(-t,v'uz) (-i)i^i|pi! n ^i^^iH'^^)- 

ZCX\Y P ■.{X\{YUZ))=[J^X, XiCP 

If Xi C X\Y, then for the trace class operator F^^^^^ and the unitary group of operator (fTS!) the 
equality is valid 

Tts+i,.. .,s+n Y\_ S\x,\i-t;Xi)Fl^x\i^)=T^^s+i,...,s+nFl'x\{X). 

Then, taking into account the equahty 

J2 (-l)|P||P|! = (-l)l^\(^^^)l 
P:(x\(yuz))=ux, 

from expression (IHOil for solution expansion (1771) we obtain 

oo ^ 

F,{t,Y) = 5^-Tr,+i,...,,+„[/i+„(-t,{y},X\y)Fi°„(X), (81) 

n=0 

where [/i+„(— t) is the (1 + n)t/i-order reduced cumulant of groups of operators ( ITSl) 
Ui^n{-t,{Y},X\Y)= (-l)l^\^^^)l%uz|(-t,FUZ). 

ZCX\Y 

Using the symmetry property of particles obeying the Maxwell-Boltzmann statistics, the equal- 
ities valid 

Tr,+i,...,,+„ Yl hiy''''^'''"%Yuzi{-t,YUZ)F;^x\iX) = 

ZCX\Y 

n s+n 

= Tr,+i,...,,+„ Y G\Y\+n-k{-t, Y,t,,..., z„_fc)Ff+„(X) = 

k=0 ii<...<j„_fe=s+l 

= Tr,+i,...,,+„ ""^.fi Yl+n-ki-t, Y,s + l,...,s + n- k)F^UX). 

Hence for the evolution operator [/i+„(— t) it holds 

f/i+„(-t, {y}, X\Y) = Y{-l)'-^^7^^as+n^k{-t, Y,s + l,...,s + n-k), (82) 
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and consequently, we derive the representation which is written down in terms of the operator 
a (an analog of the annihilation operator |T]) 

(a/)„ = Tln+lfn+l, (83) 

as the following expansion 

oo ^ n I 

nO^E^^B-D'^f^^-'Sl-Oa^F". (84) 

n=0 k=0 ^ ' 

We remark that this representation for the solution expansion for the first time is obtain in |1] 
by another method in the form 

F(t) = e''e?(-t)e-"F(0). 
Finally, in view of the validity of the equality 

^6^(-t + r)a6;(-T)F(0) = Q{-i + r) {M , ^Q{-t)F{^), 

where ■] is the commutator of operators, namely in componentwise form 

s 

( [AT, a] f){J) = J2 Tr,+i(-A/;nt(j, s + l))/.+i(r, s + 1), 
j=i 

expansion (18^ is represented in the form of the perturbation series of the quantum BBGKY 
hierarchy (!73l) 



t 

oo 



F{t) = J2 dh... / dtng{-t + h)[Ar,a]g{-ti + t2) 

n=0 { { 
. ■ ■ Qi-tn-l + tn) [U, a]g{-tn)F{Q), 



or in componentwise form 

t tn-l 
OO „ „ S 

F,{t, Y) = Y, dh... / dKi:Ts+i,...,s+nQs{-t + ti) Y^{-MU3us + 1)) X (85) 
n=0 { { n=i 

s+n—l 

xQs+li-tl + t2) • • • Qs+n-li-tn-l + Q Yl ("-^mtOn, S + n) )^,+„ (-t„) F°+„ (X) . 

in = l 

Recurrence relations (HOl) underlie of the classification of possible solution representations 
of the Cauchy problem (!73|) - (!74|) of the quantum BBGKY hierarchy. In fact, using cluster 
expansion (HOl) of the group of operators (fT5|) it is possible to construct other representations. 
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For example, solving the recurrence relations ( HOl) with respect to the cumulants of first and 
second order, we have 

2ti+„(-t,{y},x\y)= (86) 

= M-t,m,{z}) E (-i)'^'|p|! n 

ZCX\Y, P:iX\{YUZ))=[J^Xi X,CP 

Summing up the relevant terms of expression (!86|) in expansion (ITHI) similarly to the case of 
( 18T|) . we obtain the expansion representation of initial- value problem (!73|) -( [7i|) through the 
second order cumulant 

oo _ 

Fs{t,Y) = J2^^^s-,i,...,s^n Yl (-l)l''^^''"^^'2l2(-t,m,{Z})F,V(X). (87) 

n=0 ■ ZCX\Y, 

For -F(O) G 2}^{J^-u) series (1871) converges on the norm of the space 2}^{F-u) and the estimate 
holds 

II^WL,(.«)<2e1|F(0)||,.(,^, 

The statement that for initial data from the space 2}^{J^y) solution (!87|) of the Cauchy problem 
of the quantum BBGKY hierarchy (173|) - (fM|) actually follows from the equivalence of different 
representations and the existence statement of solution ( [77|) . 

4.4 The nonlinear quantum BBGKY hierarchy for marginal corre- 
lation operators 

In view of the definition of mean- value functional (130|) . for example, the dispersion of an 
additive-type observable is determined by the functional 

((A«-(A«)(t))^)(t)= (88) 

oo ^ 

= Y.-^ - {A^'Vm^Ut. 1, . . . , 1 + n) + 

n=0 



+ -T Tri,...,2+„ ai(l)ai(2)5(2+„(t, 1, . . . , 2 + n), 



n=0 



where {A^^^){t) is defined by expression (!30|) and the operators gs+n{t) are defined by expansions 
(12T|). For G £(J'-h) and ^ G ^^(J'-h) functional ([88]) exists. Following to formula §8^) we 
introduce the marginal correlation operators by the series 



oo ^ 

(t, l,...,s) = V" — Tr^+i,...,3+„ ^s+„(t, l,...,s-Fn), s > 1, (89) 



ni 

n=0 



where the operator Qs+nit, I, . . . , s + n) is defined by expansion fl2T|) over solutions f|T4|) of the 
von Neumann equations ffTTl) . According to estimate (H9|) . series (|89|) exists and the estimate 
holds: ||a(t)L.(^^) <s!(2e2)^c^E:Lo(2e')"c^ 
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Then the dispersion of an additive-type observable is defined within framework of the 
marginal correlation operators ( l89l) as follows |3] 

((A«-(A«)(t))2)(t) = 

= Tn (a?(l) - 1) + Tri,2 ai(l)ai(2)G'2(t, 1, 2). 

Thus, macroscopic characteristics of fluctuations of observable are determined by marginal 
correlation operators fl89p on the microscopic level. 

Assuming as a basis an alternative approach to the description of the state evolution within 
framework of the von Neumann hierarchy, we shall define the marginal correlation operators 
by the use of solutions of the Cauchy problem fl33|) - flM|) of the von Neumann hierarchy for 
correlation operators by formula flS^ . We note that every term of marginal correlation operator 
expansion (jHlD is determined by the [s + n)-particle correlation operator fHS]) as contrasted to 
marginal density operator expansion (172!) which is defined by the (1 + n)-particle correlation 
operator (15T!) . 

Traditionally marginal correlation operators are introduced by means of the cluster expan- 
sions of the marginal density operators fl72p governed by the quantum BBGKY hierarchy fl73l) 

P.Y = [J,X, X,CP 

where Y1p-y=\j x '-'^^^ possible partitions P of the set F = (1, . . . , s) into |P| 

nonempty mutually disjoint subsets Xj C Y. Hereupon solutions of cluster expansions fl90l) 

P.Y = \J,Xi X,CP 

are interpreted as the operators that describe correlations of many-particle systems. Thus, 
marginal correlation operators fl9T]) are cumulants (semi-invariants) of the marginal density 
operators. It is obvious that definition fl89l) follows from fl9T|) in consequence of definition fl72l) 
and relations (1521) between correlation operators of particle clusters and correlation operators 
of particles. 

The evolution of all possible states of quantum many-particle systems obeying the Maxwell- 
Boltzmann statistics with the Hamiltonian (|3]) can be described within the framework of 
marginal correlation operators governed by the nonlinear quantum BBGKY hierarchy 

j^Gsit, Y) = U{Y I G{t)) + Tr,+i ^(-A/;„,(z, s + 1)) {Gs+i{t, Y,s + l)+ (92) 

ieY 

+ ^ G\x^\{t,Xi)G\x2\{t-,X2)), 

P : {Y,s + 1) = Xi\JX2, 
ie Xi;s + le X2 

G,{t,Y)l^^ = G',{Y), s>l. (93) 

In equation fl^2]) the operator A/'(F|G'(t)) is generator of the von Neumann hierarchy fl5Bl) 
defined by 

ATiY I G{t)) = -Ms{Y)Gs{t, Y) + (94) 
+ E E $^(-A/'int(^i,^2))G|x,|(t,Xi)G|x,|(t,X2), 

P:y=XiU^2 il€Xii2eX2 
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where J2p y=Xi[JX2 possible partitions P of the set y = (1, . . . , s) into two 

nonempty mutually disjoint subsets Xi C Y and X2 C Y, and J2 p : (y, s + i) = XiU^2, 

i G Xi; s + 1 G X2 

sum over all possible partitions of the set {Y, s + 1) into two mutually disjoint subsets Xi and 
X2 such that i particle belongs to the subset Xi and s + 1 particle belongs to X2. 

We note that in case of many-body interaction potentials fl52]) hierarchy fl^^ has the form 

lGs{t,Y)=Af{Y\G{t))+ (95) 

00 n s 

+ Yl ^rs+i,...,s+i+n-ki-A^t^^\ji,---,3k,s + l,...,s + l + n-k)) X 

n=l k=l l=ji<-<jk 

^ E n G\x.\it,x.:}, 

P:(l,...,s + l + n-fc) = Ui^l^», ^»C:P 

2n(ii.---,ifc), |P| <n + i 

where we use notations accepted above. 

The nonlinear quantum BBGKY hierarchy fl92|) is derived on basis of the von Neumann 
hierarchy for correlation operators (|33|) according to definition fl89|) or on basis of the quantum 
BBGKY hierarchy for marginal density operators fl75]) according to cluster expansions fl^U]) . 
The evolution of marginal correlation operators of both finitely and infinitely many quantum 
particles is described by initial- value problem of the nonlinear quantum BBGKY hierarchy (!92|) . 
For finitely many particles the nonlinear quantum BBGKY hierarchy is equivalent to the von 
Neumann hierarchy (133!) . 

To construct a nonperturbative solution of the Cauchy problem (p2|) - (!93|) of the nonlinear 
quantum BBGKY hierarchy we first consider in capacity of initial data the marginal correlation 
operators satisfying a chaos property. In terms of marginal correlation operators fl89p a chaos 
property means that 

G°(l,...,s) = (7?(l)<5.,i, s>l, (96) 

where 6s,i is a Kronecker symbol. Taking into account the structure of a solution fl38|) of the 
von Neumann hierarchy ([3SD in case of initial data (1571) . from expansion flS^ we obtain 

00 ^ s+n 

Gs{t,l,...,s) = ^ — TTs+i,...,s+n^s+n{-t,l,---,s + n)Ylg^{i),, s > 1, 

n=0 ^' i=l 

where 2l<j+„(— t) is the (s + n)t/i-order cumulant ( l39l) of groups of operators ( ITSl) . In consequence 
of equality ( l96l) we finally derive 

00 ^ s+n 

G,(t, 1, . . . , s) = 5^ - Tr,+i,...,,+„ 2l,+,(-t, 1, . . . , s + n) J] G?(z). (97) 

n=0 ■ i=l 

Since estimate (112|) holds, series (p7|) converges provided that ||Gi(0)||£i(^) < (2e)~^. 

Thus, the cumulant structure of solution ( H^ of the von Neumann hierarchy ( 133|) induces 
the cumulant structure of expansion solution (!97|) of the initial-value problem of the nonlinear 
quantum BBGKY hierarchy for marginal correlation operators. 
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We emphasize that in case of chaos initial data solution expansion ( 177|) of the quantum 
BBGKY hierarchy ( 1731) for marginal density operators differs from solution expansion ( 197|) of 
the nonlinear quantum BBGKY hierarchy (p2!) for marginal correlation operators only by the 
order of the cumulants of the groups of operators of the von Neumann equations [28] , [H] 



oo ^ s+n 

F,{t, Y) = J2- Tr.+i,...,.+n ^i+ni-t, {Y},X\Y)l[ F^{t), s > 1, (98) 



n=0 1=1 



where the operator 2li+n(i) is (1 + n)f/i-order cumulant ( j78|l . Series ( l98l) converges under the 
condition that: ||Fi(0)||£i(-^) < e~^. 

Let us construct a nonperturbative solution expansion of the Cauchy problem (^2^- f l93p 
in the case of general initial data. According to definition (!89|) of the marginal correlation 
operators, i.e. 

G{t) = e'g{t), 

where the operator e° is defined by f lS^ and the sequence g{t) is a solution of the von Neumann 
hierarchy defined by group (148|) . i.e. 

= g{t I (7(0)), (99) 

and the equahty: g{0) = e~°G{0), we finally derive 

G(t) = e"^(t I e~°G'(0)). (100) 

Thus, a solution of the Cauchy problem of the nonlinear quantum BBGKY hierarchy for 
marginal correlation operators is defined by the following one-parametric mapping 

R3 t ^U{t \ f) = e°^(t I e'^f), 

which is defined on the space £^(J-'-h) and has the group property. On the subspaces £o(^s) C 
St^iTis), s > 1, the infinitesimal generator B{Y\f) of this group is defined by the operator 

B{Y I /) = Af{Y I /) + Tr,+i J2i-^int{t, s + 1)) {U,{t, Y,s + 1) + 

P : (y,s+l) =XiU^2, 

j e Xi ; s + 1 e X2 

where the same notations as for formula fl92|) have been used. 

To set down formula fllOOp in componentwise form, we observe that the following equality 
holds 

OO { ^ \k ^ /I 

n (e-"G(0)),,,(X.) = Y: ^Tr.,„,,...,,„,, ^ Uik-kA^ ' " " 
XiCP fc=o ■ fei=0 ^'^ 

■■■ Yl 1 \(h ^ ^' h T;G\x^\+k-kAXi,s + n+l,...,s + n + k-ki) ... 

^ A;|p|_i!(fc|p|_2 - «:|P|-i)! ' ' 



k 



|p|-i = 



• ■ ■ X G'j'x|p||+fc|p, ,(^|P|, s + n + /c - A;|p|_i + 1, . . . , s + n + /c). 
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Then in view of definitions fl83l) and (HHl), we have 



oo ^ 

Gs{t, Y) = Y,^_ Tr,+i,...,,+„ Yl {^i}' • • • ' {^|P|}) X 

n=0 ■ P:{Y,s+l,...,s+n)=[J-Xi 

X n (e-"G(0))|x.|(X,), 

X,CP 

and as a result we derive a solution expansion of the nonlinear quantum BBGKY hierarchy [I2] 

oo ^ 

= V ^ Tr,+i,„„,+„ f/i+„(t; {Y},s + l,...,s + n\ G{0)), s > 1, (101) 



n! 

n=0 



where the (1 + n)f/i-order reduced cumulant of groups ( HHj) has been introduced 

[/i+,(t; {y}, s + 1, . . . , s + n I G'(O)) = (102) 



n! 



= E(-l)'^i7;^ E 2l|p|(-t, {X,}, . . . , {X|p|}) X 

k=0 P:(l,...,s+n-fc)=UiXi 

s + n- A; + l,...,s + n-A;i)... G'pf|p||+;^.|p|_^(X|p|, s + n - A;|p|_i + 1, . . . , s + n). 

Let us indicate some properties of reduced cumulants fll02p of groups of operators psj) . In 
case of n = 0, first order reduced cumulant fll02p has the form 

m-. {y} I /) = E • • • ' n 

i.e. it is the group of operators fj48|l . Its infinitesimal generator coincides with generator of the 
von Neumann hierarchy (l36l) 



lim|(t/i(t; {Y} I /) - /,(F)) = U{Y \ /), s > 1, 

for / G £}{J^-}i) in the sense of the norm convergence of the space £}{'Hs), where the operator 
M{Y I /) is defined by formula (IMIl . 

In case of n = 1 for second order reduced cumulant (1102^ in the same sense we obtain the 
following equality 

Tr,+i lim - U2{t- {Y}, s + 1 | /) = V Ti,+^{-UUh s + 1)) (/s+i(t, Y,s + l) + 



ieY 



+ Yl f\x,\{t,X,)j]x,i{t,X,)), 



P:(y,s + l) = XiU^2, 
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where notations are used as above for hierarchy ( l92l) . and for n > 2 as a consequence of the 
fact that we consider a system of particles interacting by a two-body potential, it holds 

TTs+i,...,s+n lim ^ Ui+n{t] {Y} , s + 1, . . . ,s + n \ f) = 0. 
According to the estimate 

||2l|p|(-t, {Xi}, . . . , {X|p|})/„||^,(^^^ < |P|! elP| ||/„,||^,(^^^, 

on the space £}{1-Ls), series f llOip converges provided that max„>i IIG^ILi^t, n < (2e^)~^. 

For abstract initial-value problem for hierarchy fl92p on the space of sequences of trace-class 
operators the following statement is true \A2\. 

Theorem 6. // max„>i ||G2|Li^a/ \ < (2e^)^"'^, then fort G M a solution of the initial-value 
problem f l92p - (l93l) of the nonlinear quantum BBGKY hierarchy is determined by expansion 
f llOip . If G SyliTin) C il^(7^„) it is a strong (classical) solution and for arbitrary initial data 
Gl G H^iTin) it is a weak (generalized) solution. 

We note also that in case of many-particle systems obeying quantum statistics, i.e. many- 
particle systems of fermions or bosons the nonlinear quantum BBGKY hierarchy for marginal 
correlation operators has the form 

j^Gsit, Y) = mY\G{t)) + Tr,+i ^(-A^i.^lz, s + 1)) Y,s + 1) + 

+ 5±iG',x,l(t,Xi)G'lx,|(t,X2)), 5>1, 

P:(y,s + l) = XiU^2, 

j e Xi ; s + 1 e X2 

where M{Y\G{t)) is generator (!54l) of the von Neumann hierarchy of fermions or bosons and 
the operator Sfj^^ is defined by formula ([2]). 

We emphasize that the evolution of marginal correlation operators of both finitely and in- 
finitely many quantum particles is described by initial- value problem of the nonlinear quantum 
BBGKY hierarchy fl92|) . For finitely many particles the nonlinear quantum BBGKY hierarchy 
is equivalent to the von Neumann hierarchy ( !33l) ). 



5 The origin of the quantum kinetic evolution 

It is well known [ll|3] that in certain situations the collective behavior of many-particle sys- 
tems can be adequately described by the kinetic equations, i.e. by evolution equations for a 
one-particle marginal density operator. In this section we discuss the problem of potentialities 
inherent in the description of the evolution of states of many-particle systems in terms of a 
one-particle density operator or more exactly the problem of an equivalence of the hierarchies 
of quantum evolution equations and quantum kinetic equations. We demonstrate that in fact, 
if initial data is completely defined by a one-particle marginal density operator, then all pos- 
sible states of infinite-particle systems at arbitrary moment of time can be described within 
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the framework of a one-particle density operator and the marginal functionals of such states 
without any approximations |36], |13]. By other words the state described in terms of the 
sequence F(t) = (1, Fi(f), . . . , Fs(f), . . .) of marginal density operators can be described within 
the framework of the sequence F{t \ Fi{t)) = (1, Fi{t),F2{t \ Fi{t)), | . . .) of ex- 

plicitly defined marginal functionals Fs(t \ Fi{t)), s > 2, of the solution Fi(t) of the generalized 
quantum kinetic equation. 

Thus, in this section the origin of the microscopic description of quantum many-particle 
systems by means of kinetic equations is considered. 



5.1 The generalized quantum kinetic equation 

We consider the Cauchy problem of the quantum BBGKY hierarchy fl75]) with initial data 
F{t)\t=o = F^'^^ = {Fi{l), ■ ■ ■ ,Y[i=iFiii), ■ ■ which is intrinsic for the kinetic description 
of many-particle systems because in this case all possible states are characterized by means 
of a one-particle marginal density operator. Then we deal with initial value problem of the 
quantum BBGKY hierarchy (173|) which is not completely well-defined Cauchy problem, because 
the generic initial data F^'^\ is not independent for every unknown marginal density operator 
Fs{t,l, . . . , s), s > 1, from the hierarchy of equations. Consequently such initial- value problem 
can be naturally reformulated as the new Cauchy problem for a one-particle density operator, 
that corresponds to the independent initial one-particle density operator and the sequence of 
explicitly defined marginal functionals of the state Fs(t,l, . . . , s \ Fi{t)^ , s > 2, of the solution 
Fi{t) of such Cauchy problem [36j. 

Let us formulate the restated Cauchy problem and a sequence of the marginal functionals 
of the state which describe the evolution of all possible states of quantum many particles in an 
equivalent way as compared with the quantum BBGKY hierarchy. 

The one-particle density operator Fi{t) is governed by the following initial- value problem of 
the generalized quantum kinetic equation [36] 



^Fi(t,l) = -A/-(l)Fi(t,l)+ (103) 
at 

oo ^ n+2 

+Tr2(-Ari„t(l, 2)) ;^Tr3,...,„+2 5Ji+„(t, {1, 2}, 3, . . . , n + 2) J] Fi(t, z), 

n=0 ■ 4=1 

Fi(t,l)|,=o = i^i°(l), (104) 

where the operator A/'int(l,2) is defined by formula (^^, and the (1 + n)t/i-order generated 
evolution operator QJi+„(t), r?, > 0, is defined as follows (in case of s = 2 for F = (1, . . . , s) and 
X\Y = {s + l,...,s + n)) 

n n n-ni-...-nfc_i ^ 

5Ji+„(t,{y},x\r) = 5^(-i)'^ 5^... 7 (105) 

^ — ^ ^ — ^ ^ — ^ [n — TLi — . . . — nh)\ 

k=0 ni=l nfc=l ^ ^ 

x2li+„_„,_..._„j^(t, {y}, s + l,...,s + n-ni- ...-Uk)^ 
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k ^ s+n-ni-...-nj ^ 

|Dj I < s + 71 — ni — ■ • ■ — rij 

In expansion (11051) we denote by XId z =u Xi. sum over all possible dissections of the 

linearly ordered set Zj = {s + n — ni — ... — rij + 1, s + n — ni — ... — nj^i) on no more 
than s + n — rii — . . . — rij linearly ordered subsets and 2li+ri(^) is the (1 + n)t/i-order cumulant 

^,Ut,{Y},X\Y)= Yl (-I)I^I-^(IPI-I)! n Gmxdt,0{X.)), 

P:({r},x\y)=u,x, x,cp 

of the scattering operators fH7j) 

n 

Gait) = Gni-t, 1, . . . ,n) n > 1. 

i=l 

The series of collision integral f llOSp converges under the condition: ||Fi(t) H^ic^) < [39]. 

The marginal functionals of the state Fs{t,l, . . . , s \ Fi{t)), s > 2, are represented by the 
following expansions over products of the solution Fi{t) of the Cauchy problem fll03p - (11041) 

oo ^ s+n 

F,{t,Y\ F,{t)) = Y,^T^^s+i,...,s+n'^i+n{t,{Y},X\Y)l[F^{t,z), (106) 

n=0 ■ i=l 

where the (1 + n)t/i-order generated evolution operator QJi+„(t), n > 0, is defined by expansion 
(11051) . The marginal functionals of the state are represented by converged series (I106P under 
the condition that: \\Fi{t)\\zi[n) < 6'^^'+^^ [39]. 

We observe that the kinetic dynamics of states is described in terms of cumulants of scattering 
operators fHTI) in contrast to the evolution of states described by the quantum BBGKY hierarchy 
(!73|) . We give a few examples of the generated evolution operators QJ„, n > 1, of the lower 
orders 

Q3i(t,{y})=2ti(t,{F}), (107) 

s 

^2it, {Y},s + 1) = 2t2(t, {Y}, s + 1)- Mt, {Y}) Mt, ^,s + l). 

1=1 

It should be emphasized that in case under consideration, i.e. in case of the absence of cor- 
relations at initial time, the correlations generated by the dynamics of a system are completely 
governed by evolution operators (I105p . 

Typical properties for the kinetic description of the evolution of constructed marginal func- 
tionals of the state (11061) are induced by the properties of generated evolution operators (I105p . 
Let us indicate some intrinsic properties of the evolution operators Q3i+„ (t), n > 0, representa- 
tive for cumulants (semi-invariants) of groups of operators. 
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First of all we observe that they are solutions of the following recursive relations (the kinetic 
cluster expansions, which are in some sense analog of equilibrium virial expansions) [36] 

2ti+„(-t,{F},s + l,...,s + n) = (108) 

" n' 1 
= 7^7 '-^^i+n-mit,{Y},s + l,...,s + n-ni) V — x 

|D| < s + n — ni 



s+n—ni ^ s+n—ni 

X 



n n^2ti+|x,|(-t,Zfc,Xfc) II 2li(-t,m) 



m ^ il, . . . ,i|o| 

where ^^.^^u Xi |D|<s+n-ni ^^'^ possible dissections D of the linearly ordered 

set Z = {s + n — rii + 1, . . . , s + n) on no more than s + n — rii linearly ordered subsets. 

Since in case of a system of non-interacting particles for scattering operators fHTl) the equality 
holds: Qn{t) = I, where J is a unit operator, then we have 



2Ji+„(t) = IS, 



where Sn^ is a Kronecker symbol. Similarly, at initial time t = it is true: QJi+„(0) = /5n,o- 

The infinitesimal generator of the first-order generated evolution operator fllOSp is defined 
by the following limit in the sense of the norm convergence in the space ii^('H„) 

1 " 
hm -(2Ji(t, {!,..., n}) - /)/„ = V {-Mnt{^, j))fn, 
t-^o t ^-^ 

i<j=l 

where the operator {—J\fmt{h j)) is defined by formula f l59l) . In general case, i.e. n > 2, in 
the sense of the norm convergence on the space il^('H„) for the n-order generated evolution 
operator (llOSp it holds 

lim^QJ„(t, 1, . . .,n)fn = 0. 

Before constructing a solution of initial- value problem fll03p -( IT04|) in the space £^(7^) we 
generalize kinetic equation (I103P for particles interacting via many-body interaction potentials 
^(n)^ n > 1. In this case the generalized quantum kinetic equation has the form 

-R{t, 1) = -m)F,{t, 1) + E E 7^^^ Tr,,...,„^,(-A/;(r))(l, (109) 

n=l k=l ^ ' 

n+l 

...,k + l)QJi+„_fe(t, {1, . . . ,k + l},k + 2, . . . ,n + l)Y[Fi{t,z), 

1=1 

where QJi+„_fc(t), is the {1 + n — A;)t/i-order generated evolution operator fllOSp and notations 
dsn]), (EH]) are used. The collision integral in the generalized quantum kinetic equation f ll09p is 
defined by the convergent series under condition that ||Fi(t)||£ic^) < [36] . 
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For the sake of a comparison of the structure of various colhsion integral components in 
(11091) we give expressions of the collision integral term describing a two-body interaction and 
three particle correlations 

Tr2,3(-Ar/„?)(1, 2)2J2(t, {1, 2}, 3)F^{t, l)F,{t, 2)F,{t, 3), 

and the collision integral term describing a three-body interaction 

^Tr2,3(-<^ )(1, 2, 3)^i(t, {1, 2, 3})F^{t, l)F,{t, 2)F,{t, 3), 

where the evolution operators ^2it, {1; 2}, 3) and 5Ji(t, {1, 2, 3}) are defined by (llOSp . 

For the Cauchy problem ffTO - ffTDIj) (or ffTO - ffTuj]) ) on the space 2^{n) the following state- 
ment is true. 

Theorem 7. The global in time solution of initial-value problem (I109p - (ll04p is determined by 
the following expansion 

^) = 5Z Tr2,...,i+„ 2li+„(-t, 1, . . . , n + 1) J] F°(z), (110) 

n=0 ■ i=l 

where the cumulants 2li+„(— t), n > 0, are defined by formula ( ffSl) . // ||F°||£i(^) < (e(l + e^))^^, 
then for F° G -Co('H) a strong (classical) solution and for an arbitrary initial data F° G 
S^^lTi) it is a weak (generalized) solution. 

In section 4 the relationships of the evolution of observables and quantum states described 
in terms of marginal density operators have been considered in the general case. In case of 
initial states specified by a one-particle marginal density operator, the dual BBGKY hierarchy 
describes the dual picture of the evolution to the picture of the evolution of states governed 
by the generalized quantum kinetic equation and an infinite sequence of explicitly defined 
functional of a solution of such evolution equation. In fact, the following equality is true 

(5(t),FW) = (S(0),F(t|Fi(t))), (111) 

where the initial state F'^ is given as above, the sequence B{t) defined by expansion ( l63i) and 
F{t I Fi{t)) = (1, Fi{t), F2(t I Fi{t)), . . . , Fs{t I Fi{t)), . . .) is a sequence of marginal functional 
of the state (11060 with the first element Fi{t) given by series (IllOp . 
To verify equality (II lip we transform functional {B{t)\F^^^) as follows 

oo ^ 

{B{t), fW) = J2- Tri,...,s i?°(l, . . . , 5) X (112) 

oo ^ s 

2li+„(-t, {Y}, s + l,...,s + n)l[ F°(z), 

n=0 ^' i=l 

where the (1 -|- n)t/i-order cumulant 2li+„(— t, {Y}, s + 1, . . . , s -|- n) is defined by (ITHI) . For F^ G 
2^(71) and G -C('Hs) obtained functional (I112p exists under the condition ||F°||£i(-^) < e^^. 
Then we expand the cumulants 2ti+„(— t) over the new evolution operators QJi_|_„(t), n > 0, 
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into the kinetic cluster expansions f llOSp . Representing series over the summation index n and 
the sum over the summation index rii in functional (11121) as the two-fold series and identifying 
the series over the summation index ni with the products of one-particle density operators 



oo s+n 

^ Tr,+„+i,...,,+„+„, ^ XI n ^ '^^^^^ 

"1=0 D:Z = Ufc^fe, n<---<jlD|=l ^fcCD I ^1' 

D| < s + n 



s+n n+s+n-i s+n 

x2ti+ix,i(-t,zfc,Xfc) n 2ti(-t,o n ^i°(j)=n^i(^'')' 

I = 1, i=i i=l 

I il,..., ipi 

we transform functional flll2p to the form in terms of marginal functionals of the state fll06p . 
Thus, equality fillip holds. 

In a particular case of the additive-type marginal observables B^^\t) given by ( 166|) equality 
f illip is reduced to the form 

(S«(t),F(^)) = Triai(l)Fi(t,l), 
where the one-particle marginal density operator Fi{t) is represented by the expansion 

oo ^ n+1 

^1^^' ^^ = T.^_ Tr2,...,i+„ ^i+n{-t, 1, . . . , n + 1) n F{'(z), 

n=0 ■ 1=1 

which coincides with solution f lllOp of the Cauchy problem f ll03p -( |T0^ . Hence for additive- 
type marginal observables the generalized quantum kinetic equation fllOSp is dual to the dual 
quantum BBGKY hierarchy fl57j) with respect to bilinear form flTTj) . 

Let us derive the evolution equation, which satisfies obtained expansion for the one-particle 
marginal density operator Fi{t). Taking into account equality flT3|) and observing the validity 
of equahty fl7^ for cumulants of groups ffT^ . we differentiate over the time variable in the sense 
of pointwise convergence on the space £^(7/) obtained expansion for Fi{t). As result it holds 

|Fi(t,l) = -M(l)Fi(t,l)+ (114) 



oo ^ n+2 

+Tr2(-Mnt(l, 2)) ^ -Tr3,...,n+2 2li+„(t, {1, 2}, 3, . . . , n + 2) J] Fff^). 



n=0 1=1 



In the second summand in the right-hand side of equality f lll4p we expand cumulants (1781) 
of groups ( JT5l) into kinetic cluster expansions (llOSp and represent series over the summation 
index n and the sum over the summation index rii as the two-fold series. Then the following 
equalities take place 



n+2 

0/ 



J2 ^ Tr2,...,„+2(-A/;nt(l, 2))2ti+„(t, {1, 2}, 3, . . . , n + 2) J] F{ 

n=0 ■ i=l 
oo ^ n I 

= 5" ^ Tr2,...,„+2(-Ari„t(l, 2)) V gJi+„_„,(t, {1, 2}, 3, . . . , n + 2 - m) X 

n=0 ni=0 ^ '-^ 
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n+2—ni 2+n— ni ra+2 

D:Z=U,Xi ' n7^...7^i|D|=l ^;CD I 'I' m = 1, i=l 

m ^ il, . . . 



Tr2(-Ari,t(l, 2)) 5^ - Tr3,...,n+22Ji+.(t, {1, 2}, 3, . . . , n + 2) — x 



n=U ■ ni=0 D:Z'=[JiXi 

n+2 n+2 n+2+ni 

^ E n n n 

where Z = (n + 3 — ni , . . . , n + 2) and Z' = (n + 3, . . . , n + 2 + ni) are hnearly ordered sets and 
the notations accepted above are used. Consequently, in case of s = 2 applying to the obtained 
expression formula f lll3p . from equality f lll4p we derive 

|Fi(t,l) = -M(l)Fi(t,l) + 



oo ^ n+2 

+Tr2(-A/'int(l, 2)) - Tr3,...,n+2 5Ji+„(t, {1, 2}, 3, . . . , n + 2) J] Fi(t z) 



n! 

71=0 i=l 

Constructed identity for a one-particle (marginal) density operator defined by expansion fill 01) 
we treat as the evolution equation, which governs the one-particle states of many-particle 
quantum systems. 

We remark that one more approach to the derivation of the generalized quantum kinetic 
equation consists in its construction on the basis of dynamics of correlations [H], |28j . 

Thus, in case of initial data F'^'^^ = (F{'(1), . . . , ni=i-^i°(^)i • • ■)? which is completely charac- 
terized by the one-particle marginal density operator F^, solution ( 1631) of the Cauchy problem 
(!57|) - (l58|) of the dual quantum BBGKY hierarchy for marginal observables and a solution of 
the Cauchy problem of the generalized kinetic equation fll03l) - fll04p together with marginal 
functional of the state ( I106P give two equivalent approaches to the description of the evolution 
of quantum many-particle systems. 

We note also that in case of many-particle systems obeying the Fermi-Dirac and Bose- 
Einstein statistics the generalized quantum kinetic equation has the form [37] 

^Fi(t,l) = -Ar(l)Fi(t,l) + 

oo ^ n+2 

+Tr2(-Ari„t(l, 2)) ;^Tr3,...,n+2 2Ii+n(t, {1, 2}, 3, . . . , n + 2)5±+2 H 

n=0 ' i=l 

where the operator Sf is defined by (l2l), and the marginal functional of the states are repre- 
sented by the series 

oo _ s+n 
n=0 ■ i=l 

where the (1 -|- n)th-oidei generated evolution operator QJi+„(if:), ?7, > 0, is determined by 
expansion (llOSp . 
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5.2 The marginal functionals of correlations 

Within the framework of the description of states in terms of marginal correlation operators ( 19T|) 
the marginal correlation functionals Gs{t,Y \ Fi(t)), s >2, are represented by the expansions 
similar to f ll06p . namely 

Gs{t,Y I = (115) 

CXD ^ S+n 

= J2- Tr.+i,...,s+n2Ji+„(t, 0{{Y}),s + l,...,s + n)Y[ F^{t, z), 



ni 

n=0 i=l 



where it is used the notion of declasterization mapping fl43|) . Hence in contrast to expansion 
(11061) the n term of expansions (IllSp of marginal correlation functionals Gs{t,Y \ Fi{t)) is 
governed by (1 + n)t/i-order evolution operator ( llOSp of the (s + n)t/i-order cumulants of the 
scattering operators, for example, as compared with (I107p the lower order evolution operators 
QJi+,i(t, 9{{Y}), s + 1, . . . , s + n), n > 0, have the form 

^,{t,9i{Y}))=^si-t,9{{Y}), 
^2it,ei{Y}),s + l) = 

s 

= 2i,+i(-t, e{{Y}), s + i)- 2t,(-t, e{{Y})) J2 M-t, ^,s + l), 

i=l 

and in case of s = 2, it holds 

^^{t,9{{l,2})) = g2{t,l,2)-I, 

where ^2(^5 1,2) is scattering operator fHTl) . 

We indicate that expansions (11061) of marginal functionals of the state are nonequilibrium 
analog of the Mayer-Ursell expansions over powers of the density of equilibrium marginal density 
operators. 

Within the framework of the description of states by marginal functionals of the state (11060 
the average value, for example, of the additive-type marginal observable B^^^ = (0, ai(l), 0, . . .) 
is given by the functional 

(5W>(t)=Triai(l)Fi(t,l), (116) 

i.e. it is defined by a solution of the generalized quantum kinetic equation (I103p . or in general 
case of the s-ary marginal observable i?*^*^ = (0, . . . , 0, as{l, . . . , s), 0, . . .) by the functional 

= i. Tri_, a,(l, . . . , s)Fsit, 1, . . . , s | F^{t)), s > 2, 

where Fs{t,l, . . . , s \ Fi{t)) is the marg inal functional of the state ffT06D . For B^'^ G £(Jk) and 
Fi{t) G il^('H) these functionals exist. 

The dispersion of an additive-type observable is defined by a solution of the generalized 
quantum kinetic equation (11031) and marginal correlation functionals (IllSp as follows 



;(i?w-(s«>(t))0(t) = 

= Tn (a?(l) - (S«>'(t))Fi(t,l) + Tri,2ai(l)ai(2)G2(t,l,2 | Fi(t)), 
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where the functional {B^^^)(t) is determined by expression f lll6p . 

We note that the dispersion of observables is minimal for states characterized by marginal 
correlation functionals flllSp equals to zero, i.e. from macroscopic point of view the evolution 
of many-particle states with the minimal dispersion is the Markovian kinetic evolution. In 
fact functionals (11151) or (11061) characterize the correlations of states of quantum many-particle 
systems. We illustrate close links of functionals (IllSp and (I106p in the following way 

F2(t, 1, 2 I Flit)) = Flit, l)Fi(t, 2) + G2it, 1, 2 | 

Basically this equality gives the classification of all possible currently in use scaling limits. In 
the scaling limits it is assumed that chaos property of initial state preserves in time, i.e. the 
scaling limit means such limit of dimensionless parameters of a system in which the marginal 
correlation functional G2it, 1,2 | vanishes. According to definition (IllSp . it is possible, 

if every finite particle cluster moves without collisions. 

Finally we point out the relationship of generalized quantum kinetic equation (11031) and the 
specific quantum kinetic equations. The last can be derived from the generalized quantum 
kinetic equation in the appropriate scaling limits [22] or as a result of certain approximations. 
Let us consider first two terms of expansion (11060 . If an interaction potential in ([3]) is a bounded 
operator and fs+i G Q^iTis+i), then for the second-order cumulant of scattering operators (jUj) 
an analog of the Duhamel equation holds 

2l2(t,{y},s + l)/,+i = / drgsi-r,Y)gii-r,s + l)J2i-MntituS + l))x 

n=i 

s+l 

x^,+i(t - t, y, s + 1) JJ e;i(r, Z2)/s+i, 

42 = 1 

and, consequently, for the second-order evolution operator 2J2(^, s + 1) we have 



QJ2(t, {Y}, s + l)/,+i = (2t2(t, {y}, s + l)- 2ti(t, {Y}) J2 Mt, ^us + l))/,+i 

ft s 

= / dTgsi-T,Y)gii-T,s + i){y2i-Afintiti,s + i))gs+iir-t,Y,s + i)- 

"^0 n=l 

s s+l 

-g,ir - t, Y) ^(-ATi^tl^i, s + l))^2(r - t, zi, s + 1)) J] g^ir, Z2)/.+i. 



«1 = 1 42 = 1 

Observing that in the kinetic (macroscopic) scale of the variation of variables ^2j the groups of 
operators ([7]) of finitely many particles depend on microscopic time variable e~^t, where e >0 
is a scale parameter, the dimensionless marginal functionals of the state are represented in the 
form: Fsi^e~^t,Y \ Fi(t)). Then in the limit £ — )■ the first two terms of the dimensionless 
marginal functional expansions (11060 

s 

gsie-H,Y)l[Fiit,t)+ 

i=l 
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+ / dTg,i-T,Y)TTs+i{y2i-^UH,s + l))gs+ii6-H,Y,s+l)- 
"^0 n=l 

s s+1 

-Gsie-h, Y) Y,i-^^nt{^l,s + l))^2(£"'t, ^l, s + 1)) H ^2)^1 (i, ^2) 

ii=l 42=1 

coincide with corresponding terms constructed by the perturbation method with the use of the 
weakening of correlation condition by Bogolyubov [1]. Thus, in the kinetic scale the collision 
integral of the generalized kinetic equation fllOSp takes the form of Bogolyubov's collision in- 
tegral PP and we observe that in a space homogeneous case the collision integral of the first 
approximation has a more general form than the quantum Boltzmann collision integral. 



5.3 On quantum kinetic equations in case of correlated initial data 

We have proved that in case of initial data which is completely defined by a one-particle density 
operator, all possible states of infinite-particle systems of bosons or fermions at an arbitrary 
moment of time can be described within the framework of a one-particle density operator and 
explicitly defined functionals of such operator without any approximations. One of the advan- 
tages of this approach is the possibility to construct the kinetic equations in scaling limits in 
the presence of correlations of particle states at initial time, for instance, correlations char- 
acterizing the condensed states of interacting particles obeying Fermi-Dirac or Bose-Einstein 
statistics [E]. 

Let us consider initial data 

F(t)|,=o = {F?{l),h{l, 2)^0(1)^0(2), . . . , /i„(l, n)n" , • • • ) ' 

■*--*-i=i 

where the operators hn G Il^('Hn), n > 2, are specified initial correlations. Such initial data 
is typical for the condensed states of quantum gases, for example, the equilibrium state of the 
Bose condensate satisfies the weakening of correlation condition with the correlations which 
characterize the condensed state [3]. 

In the case under consideration the kinetic cluster expansions of cumulants of groups, i.e. 
recurrence relations (11081) . take the form 

s+n 

2li+„(-t, {Y}, s + l,...,s + n)/ii+„({y}, s + l,...,s + n)Y[Qii{t,i) = (117) 

1=1 

n I s+n—ni 

= Y. f ., ^i+n-nAtAY},s + l,...,s + n-ni) V V 

^ — ^ n — Til ! ^-^ ^-^ 

|Db+„| ^ |De+n| 

n ^2ti+|x,|(t,^fc,Xfc) J] /ii+|x,|(zfc,Xfc)2li(t,4) J]2li(t,j), 

fc=i ' k=i jez 

where the notations of formula (11081) are used. In terms of new evolution operators (the 
(1 + n)th-oTdeT scattering cumulants) 

s+n 

2ti+„(t, {Y}, X\Y) = 2li+„(-t, {Y}, X\y)/ii+„({F}, X\Y) J] 2li(f, t), 

i=l 
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the solutions S5i+„(t, {!"}, n > 0, of recurrence relations (11171) are represented by the 

expansions 

ri n n~ni-...-nk-i 

<s.+,.(un.Y\n = "!E(-i)'E- E („_„,_ ("s) 

fc=0 ni = l nfe=l ^1 

x2li+„_„^_..._„^(t, {y}, s + 1, . . . , s + n - rii - . . . - Hk) X 

Dj I < s + n — ^^l — ■■■ — no- 
where J2Dj-Zj=[Ji , Xi. possible dissections of the hnearly ordered set Zj = 

[s + n — rii — . . . — rij + 1, . . . , s + n — rii — . . . — rij^i) on no more than s + n — rii — . . . — rij 
linearly ordered subsets. For example, 

s 

&,{t, {Y}) = 2ti(t, {Y}) = J]2ti(t, {). 

Therefore the one-particle density operator Fi{t) is governed by the following generalized 
quantum kinetic equation 

|Fi(t,l) = -A/-(l)Fi(t,l)+ (119) 

oo ^ n+2 

+Tr2(-Mnt(l, 2)) ;^Tr3,...,n+2 &i+n{t, { 1 , 2} , 3, . . . , n + 2) J] Fi(t, ^), 

n=0 ■ 1=1 

where the operator A/'int(l,2) is defined by formula ( ffSl) . and the (1 + n)t/i-order generated 
evolution operator 0i+„(t), n > 0, is defined by expansion (IllSp . Correspondingly the marginal 
functionals of the state are represented by the following expansions 

oo s+n 

F,{t, Y I Fi(t)) = E Tr,+i,...,,+„ 0i+„(t, {Y}, X \ F) J] F,{t, t). (120) 

n=0 ■ i=l 

Thus, the coefficients of generalized quantum kinetic equation (11191) and generated evolu- 
tion operators (11181) of marginal functionals (11201) are determined by the operators of initial 
correlations. 



6 On scaling limits of hierarchy solutions 

The current point of view on the problem of the derivation of quantum kinetic equations from 
underlaying many-particle dynamics consists in the following. Since the evolution of states 
of infinitely many quantum particles is generally described by a sequence of marginal density 
operators governed by the quantum BBGKY hierarchy, then the evolution of states can be 
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effectively described by a one-particle density operator governed by the quantum kinetic equa- 
tion in a suitable scaling limit [221123] . In this section the mean field (self-consistent field) 
scaling asymptotic behavior of stated above solutions is established. The constructed asymp- 
totics are governed by the quantum Vlasov hierarchy for limit marginal density operators and 
limit solution of the generalized quantum kinetic equation [M], the nonlinear quantum Vlasov 
hierarchy for limit marginal correlation operators and the dual quantum Vlasov hierarchy for 
limit marginal observables, respectively [13]. In case of initial data satisfying a chaos prop- 
erty [1], which means the absence of correlations at initial time, the constructed asymptotics 
are governed by the Vlasov quantum kinetic equation. Moreover, within the framework of 
the description of the evolution by the dual quantum BBGKY hierarchy or by the generalized 
quantum kinetic equation it is possible to construct the kinetic equations in scaling limits in 
case of presence of correlations at initial time, for instance, correlations which characterize the 
condensate states [38] . 

6.1 A mean field limit of a solution of the dual quantum BBGKY 
hierarchy 

We consider the n-particle system with the Hamiltonian 

n n 

H^ = J2Ki^+e $(21,^2), (121) 

i=l n<i2 = l 

where e > is a scaling parameter. At first we construct the mean field scaling limit of a 
solution of initial-value problem (I57p -f l58p of the dual quantum BBGKY hierarchy. 
Let for initial data 5° e ^{'Hs) there exists the scaling limit 6° G £(7/5), i.e. 

w*- lim(e-^S° - 6°) = 0, (122) 

then for arbitrary finite time interval there exists the mean field limit of solution fl63|) of the 
Cauchy problem fl57j) - fl58l) of the dual quantum BBGKY hierarchy in the sense of the *-weak 
convergence of the space S.{l-is) 

w*- Wmie-'BJt) - bJt)) = 0, s > 1, (123) 

e-i-O 

which is defined by the expansion 

s—l * s 
bsit,Y) = J2 fdti... I dtng^,{t-h) J2 A/'int(n,Jl)e^°-l(tl-t2) (124) 

n = g g n^jl=l 

s 

••• Qs-n+li'tn-1 — tn) M\vA.{im jn) ^ 

"^71 7^ jn — li 
in,jn 7^ (il, . . ■ , in-l) 

X^°_„(t„)Cn(n(jl,---,Jn)), 
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where the following notation of the group of operators of noninteracting particles is used 

Gs-n+li^n-l - tn) = Q^s-n+l{tn-l " ^n, ^\(jl, • • ■ , jn-l)) = 
= IJ Glitn-l-tnJ)- 

jey\(ji,...,jn-i) 

Before to prove this statement we give some comments. If 6(0) G the sequence 

b{t) = (60, bi{t), . . . , bs{t), . . .) of limit marginal observables fll24p is a generalized global solution 
of the initial-value problem of the dual quantum Vlasov hierarchy 

-b,{t,Y) = J2^i^)bs{tX)+ ^^nt{ji,32)bs-i{t,Y\ij,)), (125) 

bs{t) \t=o= b',, s>l. (126) 

This fact is proved similar to the case of an iteration series of the dual quantum BBGKY 
hierarchy [iO]. It should be noted that equations set (11251) has the structure of recurrence 
evolution equations. We make a few examples of the dual quantum Vlasov hierarchy ( 11250 in 
terms of operator kernels of the limit marginal observables 



a . . 1 
dt 



i T^hit, qi; q'l) = --(-Ag, + A,;)6i(t, qi, q[), 



d 1 ^ 

i ^^2(t, gi, q2] q'l, ^2) = "2 + ^g^)^2(t, ^i, ^2; q'l, q'2) + 

i=l 

+ mq'i - ^2) - "^(Qi - Q2)) {hit, qu q[) + biit, q2; q'2)) . 

Let us consider a particular case of observables, namely the mean field limit of the additive- 
type marginal observables. In this case solution ( l63ll of the dual quantum BBGKY hierarchy 
( 1^ has the form (IMj) . If for the additive-type observables B^^\0) = (0, ai(l), 0, . . .) condition 
(I122p is satisfied, i.e. it holds 

w*-lim(e-^ai(l) -fe;(l)) = 0, 
then according to statement (I123p . we have 



w*-lim(6-S«(t)-6«(t))=0, 



where the limit operator b^^\t) is defined by the expression 

t ts-2 



b^^\t,Y) = Jdh... J dt^^iG'lit-h) J2 A/'int(^i,Ji)6^°-i(ti-t2) (127) 

n7^ii=i 

s 

■■■ Q2{^s-2-ts-l) ^ A/'int(^s-l,is-l) X 

is-l js-i = 1, 

is-l,js-l ^ (jl, ■ ■ ■,js-2) 

Xg',its-l)b'liY\Uu...,Js-l)), S>1, 
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as a special case of expansion (11241) . We give examples of expressions f ll27p 
6«(t,l) = ^i(t,l)6?(l), 



t 



b', 'it, 1,2) 



[drflg^it- T, i) NU^, 2) ^Ar. j) b%j). 
i i=i 



To establish the relationship of constructed mean field asymptotic behavior of marginal 
observables with asymptotic behavior of marginal states we consider initial data satisfying the 
factorization property or a chaos property [1] , which means the absence of correlations at initial 
time. For a system of identical particles, obeying the Maxwell-Boltzmann statistics, we have 

s 

F(t)u=F(^)=(F{'(i),...,n^i°«>---)- 

i=l 

This assumption about initial data is intrinsic for the kinetic description of a gas, because in 
this case all possible states are characterized only by a one-particle marginal density operator. 
Let 

hm II eFO-jO 11^,(^^ = 0, (128) 



then the limit of initial state satisfies a chaos property too 

s 

/('=)^(/°(i),...,n/i°w,---)- (129) 

1=1 

For b{t) G 2^{J^y) and G £,^(71), under the condition < 7, the mean field limit of 

mean value functional ( I7T]) exists and it is determined by the expansion 



CXD S 

mj-'^'^) = Y.- tt,,...,s bs{t, 1, ...,.) n fii^- 



si 

s=0 i=l 



In consequence of the validity of the following equality for the limit additive-type marginal 
observables (I127P 

oo ^ s 

{b^'\t), /W) = 5^ - Tn,...,. 6«(t, 1, . . . , 5) n = (130) 

s=0 ■ 1=1 

= Tri6?(l)/i(t,l), 

where the operator b^f^\t) is given by expansion (I127P and the one-particle limiting density 
operator fi(t, 1) is determined by the series 

* 1 

f,{t,l) = J2 /^^l--- / dtnTT2,...,n+lY[gi{-t + h,H)x (131) 
n=0 Q Q ii=l 

2 n 

x{-Arintii,2))Y[gii-h + t2,3i)...Y[gi{-tn + tn,in) X 

jl = l in = l 

n n+1 n+1 



Theory of quantum evolution equations 



49 



we establish that operator f ll3ip is a solution of the initial-value problem of the Vlasov quantum 
kinetic equation 

1) = -Afil)m 1) + Tr^i-AfUl, 2)) 2), (132) 

m\t=o = f?. (133) 

Correspondingly, a chaos property in the Heisenberg picture of the evolution of quantum 
many-particle systems is fulfilled. This fact follows from the equality for the limit k-aiy marginal 
observables, i.e. b^''\0) = (0, . . . , bl{l, . . . , A;), 0, . . .), 

oo _ s 

{b^'\t), /(^)) = 5^ - Tn,...,. ftf (t, 1, ...,.) n m^) = (134) 

s=0 ' i=l 

k 

= ;^Tri,...,,6°(l,...,fc)f]/i(t,^), k>2, 

i=l 

where the limit one-particle marginal density operator /i(t, i) is defined by expansion (11311) and 
therefore it is governed by the Cauchy problem (ll32p -( |T33l) . 

Thus, in the mean field scaling limit an equivalent approach to the description of the kinetic 
evolution of quantum many-particle systems in terms of the Cauchy problem (I132p - (I133I) of the 
Vlasov kinetic equation is given by the Cauchy problem (I125p - (ll26p of the dual quantum Vlasov 
hierarchy for the additive-type marginal observables. In case of the A;-ary marginal observables 
a solution of the dual quantum Vlasov hierarchy ( I125P is equivalent to the preservation of a 
chaos property for /c-particle marginal density operators in the sense of equality (11340 . 

6.2 A mean field limit of a solution of the quantum BBGKY hier- 
archy 

Within the framework of the description of the evolution in terms of states the scaling mean field 
limit of a nonperturbative solution of the Cauchy problem (173!) - fl74l) of the quantum BBGKY 
hierarchy in case of arbitrary initial data is stated as follows. 
Let for initial data Ff, /° G 2}{1-Ls) it holds 

lim ||e^Ff (1, . . . , s) - /°(1, . . . , s) = 0, (135) 

then for any finite time interval for solution fl77|) there exists the mean field limit 

lim \\e'F,{t, 1, . . . , s) - /,(t, 1, . . . , s) ||^,(^^) = 0, (136) 

where the limiting operator fs{t) G 2}{1-Ls) is determined by the series 

„ * tn-l 

OO „ „ S 

s+l,...,s+n 

n=0 { { ji=l 
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^(-A/;nt(^l, S + 1)) n + ^l) • • • 



s+1 

X 

ii = l h = l 

s+n—1 s+n—1 s+n 

JJ ^l(-^n-l + tnJn) (~-^int(^n, S + n)) JJ Qfs+n^ 

which is norm convergent on the space on finite time interval. Limiting marginal 

operators (11371) are governed by the limiting BBGKY hierarchy known as the quantum Vlasov 
hierarchy 

, s s 

j^fsit) = J2^-m^))fs{t) + ^ Tr,+i(-A/;„t(2, s + s > 1. (138) 

1=1 i=l 

The validity of this statement is the consequence of the following formulas on an asymp- 
totic perturbation of cumulants of groups. If fs G £^(^5), then for the first-order cumulant 
2ti(-t, {Y}) the equality holds 

s 

lim||(2li(-t,{r})-n^iH,j))/sLi(«,) = 0, (139) 

and in general case for n > 1 we have 

hni II (11 2li+„(-t, {Y}, s + 1,..., s + n)- (140) 

* s s s+1 

- dh... / dt^Y[gi{-t + h,ji)J2i-^ini{H,s + i))Y[gi{-h + t2,h)... 

•'^^-'^ '1^-^ 

s+n—1 s+n—1 s+n 

... Y\_ Ql{-tn-l + tn, in) ("A/'int (^n , ■§ + n) ) TT ^1 (~^» ^ ^^)) /s+" 1 1 d m , , ^ = 0" 

for finite time interval. 

Let us consider initial data satisfying the factorization property or a chaos property (I129p . 
which means the absence of correlations at initial time. We observe that, if the initial data 
satisfy the chaos property for particles obeying the Maxwell-Boltzmann statistics, then solution 
expansion (I137p of initial- value problem (I138p - (ll29p of the Vlasov hierarchy is possessed of the 
same property 

s 

Ut,l,...,s) = \{h{t,3), s>2, (141) 

where fi{t,j) is defined by series (11370 in case of s = 1 and initial data (I129p . i.e. (I13ip . 
which for bounded interaction potential is the norm convergent on the space £^(7{) under the 
condition 

t<to = {2\ms,in,)\\fihHn)r'- 
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In other words, under the condition on initial data G ii^('Hs) 

s 

lim||6^Ff(l,...,5)-n/i°(j)Li(«,)=0, 
on finite time interval for solution (I77|l of the quantum BBGKY hierarchy it holds 



lim||6^F.(t,l,...,.)-n/i(t,j)L(^.) = 0, 

where the one-particle limiting density operator fi{t, 1) is determined by series (11311) which is 
a solution of the Cauchy problem fll32p - fll33p of the Vlasov quantum kinetic equation. 

In case of pure states, i.e. if the limit one-particle density operator /i(t) = li/jf) {i/jt\ is a 
one-dimensional projector onto a unit vector li^t), we have 

lim\\e^F,{t)-\ijt){^tr\\s,.^^^^ = 0, 

where \ifjt) is a solution of the Hartree equation, which in terms of the kernel fi{t,q,q') = 
ipit, q)i^{t, q') of the marginal operator fi{t) has the form 

i-ij{t,q) = --Agi:{t,q) + J dq'^q - q')\ij{t, q')\^ij{t, q) . 

If the kernel of the interaction potential is the Dirac measure = S{q), then we derive the 
cubic nonlinear Schrodinger equation 

= -^A,^(t,g) + \t/j{t,q)\^^ljit,q). 
We observe that in the general case of many-body potentials with the scaled Hamiltonian 

n n n 

1=1 k=2 ii<...<ik=l 

the Hartree equation takes on form 
d 1 

oo _ „ n+1 

+ dq2...dqn+l^^''+'\ql,...,qn+l)l[mt,q^)\'^lJ{t,ql). 

n=l ■ 1=2 

We note also that in case of many-particle systems obeying quantum statistics statement 
f ll36p is true for operators from the corresponding spaces £,^{1-Lf). For initial data satisfying a 
chaos property 

/w^(/{'(i),...,5±n' 

where the operator Sf" is defined by ([2]), the Vlasov quantum kinetic equation has the form 

1) = -Ar(l)/i(t, 1) + Tr2(-A/;nt(l, 2))52^/i(t, l)/i(f, 2), 
and consequently for pure states of fermions we derive the Hartree-Fock equation. 
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6.3 On a mean field asymptotics of dynamics of correlations 

We give some comments on the existence of mean field scaling limits of the constructed solution. 
The mean field limit gs{t, 1, s), s > 1, of marginal correlation operators fllOip exists 



provided that 



lim||e'G,(t) - c/,(t)||^,(^^^ = 0, s > 1, 



lim||e^G°-(7°||^,(^^) =0, s>l, (143) 



and it is governed by the nonlinear Vlasov quantum hierarchy 

%{t, Y) = J2(-X(i))g^(t, Y) + Tr,+i ^(-A/;,t(^, s + 1)) X (144) 



ieY ieY 



X 



{gs+i{t,Y,s + l)+ Yl gix,\{t,Xi)g\x,\{t,X2))., 



P: (Y,s + 1) = XiU^2, 
ie Xi;s + l£ X2 

where we use the same notations as for hierarhy ( l92l) . 
If initial data satisfies chaos property, then we establish 

lim||e^a(t)||^,(^^) = 0, s>2, (145) 

since solution expansions ( 1971) for marginal correlation operators are defined by the (s + n)th- 
order cumulants as contrasted to solution expansions (177|1 for marginal density operators which 
defined by the (1 + n)t/i-order cumulants and in the same way as statement fll40p the equality 
holds 

\\\^S+n{-t, 1, . . . ,S + n)f,+J\al(n, ^ -> = 0. 

In the case of s = 1 provided that fll43p we have 

lim||eG'i(t) - 9i{t)\\^^(^^^ = 0, 

where for finite time interval the limiting one-particle marginal correlation operator gi{t, 1) is 
given by the norm convergent on the space £^('H) series 

gi{tA)= (146) 

CO p p z 

= J2 dh--- / dtnTT2,...,n+lQl{-t + h,l){-Mnt{l,2))Y[gi{-ti + t2,Jl)... 
n=0 { { n=l 

n n n+1 n+1 

... n Gl{-in+tn,ln) (-A^i^t ( A;„, n + 1) ) W Gl{-tn,3n)\{gi{T). 

i„=l kn = l jn = l «=1 

which obviously coincides with iteration series f ll3ip of the Vlasov quantum kinetic equation. 
In view of the validity of limit fll45p from the Vlasov nonlinear quantum hierarchy fll44p we also 
conclude that limit one-particle marginal correlation operator (11460 is governed by the Cauchy 
problem of the Vlasov quantum kinetic equation (I132p - (ll33p . 

Therefore the Vlasov nonlinear quantum hierarchy (I144p describes the evolution of initial 
correlations. 
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6.4 On scaling limits of the generalized quantum kinetic equation 

In this section we consider the relationship of the specific quantum kinetic equations with 
the generalized quantum kinetic equation. First we construct the mean field (self-consistent 
field) asymptotics of solution (11 101) of initial-value problem of the generalized quantum kinetic 
equation for a system with the Hamiltonian (11211) . i.e. 

lF,{t,l) = -Af{l)F,{t,l)+ (147) 

oo ^ n+2 

-F5Tr2(-A/;„t(l, 2)) ^ ^Tr3,...,n+2 2Ji+„(t, {1, 2}, 3, . . . , n + 2) J] Fi(t, i), 

n=0 ' i=l 

F,{t,l)\t=o = F^{l), (148) 

and marginal correlation functionals (I115p . 

If there exists the mean field limit /° G ii^('H) of initial data (I148P in the following sense 

!i-|h^°-/l'IU(«) = 0' 

then for arbitrary finite time interval there exists the limit fi{t) of solution (IllOp of the gener- 
alized quantum kinetic equation (11471) 

lim||eFi(t)-A(t)||^,^^^^=0, (149) 

where the one-particle limiting density operator fi{t, 1) is determined by series (11311) which is 
a strong solution of the Cauchy problem of the Vlasov quantum kinetic equation 

^/i(t, 1) = -Ar(l)/i(t, 1) + TT2i-Mnt)il,2)Mt, l)/i(t, 2), 
fi{t)\t=o = f?. 

Taking into account equality (I149p . for marginal functionals of the state (11060 we establish 

s 

Jim \\e'Fs{t,Y\ F,{t)) ~l[fi{t,j)\\^,^^^) = 0, 

j=i 

where fi{t) is defined by series (11311) . This equality means that in the mean field scaling limit 
initial chaos property preserves in time. 

The validity of this limit statement is the consequence of formulas (11391) and (11401) on an 
asymptotic perturbation of cumulants of groups and definition (I105p of the generated evolution 
operators Q3i+„(t), > 0, of marginal functionals of the state (I106p . Indeed, if fs G £^(7/^), 
then the equality is correct 

;im||(^.(t,F)-/)/.||,,(^^) = 0. 

Further for the first-order generated evolution operator 2Ji(t, {Y}) the following equality holds 

hm||(2J,(Uy})-/)/.^^^^=0, 
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for any finite time interval, and in the general case n > 1, we have 

lim ||^2J,+„(t, {r},x \ Y)u4^,^^^^^^^ = 0. 

According to this formula on an asymptotic perturbation of evolution operators (11051) . we 
establish the mean field asymptotics of marginal correlation functionals (11151) 

lim \\e'Gsit, Y \ F,{t))\\^,^^^^ = 0, s > 2. 

In case of initial states involving correlations for evolution operators (IllSp in the mean field 
limit the equality holds 

lim ||(Si+„(t, {¥}, s + 1, . . . , s + n)fs+n\\^i,.y ^ = 0, n>l, 
and in case of the first-order generated evolution operator (IllSp we have respectively 

s s 

iim||(0i(t,{y}) - n^iH,^i)^i(m) n^i(^'^2))/.L.(^^) =0. 

11 = 1 12 = 1 

Thus, the mean field asymptotics of marginal functionals ( I115p has the form 

s s s 

lim||e^F,(t,y I F,{t)) - l[g,{-t,t,)h,{{Y}) l[0l{t,^2)l[fl{t,J)\\^,^^^^=O, 

il=l 12=1 i=i 

that means the propagation of initial correlations in time in the mean field limit, and the limit 
one-particle density operator satisfies the modified Vlasov quantum kinetic equation 

±f,(t,i) = -m)fiit,i)+ (150) 

2 2 

+^T2{-Mnt){l,2)l[g^{-t,^^)hi{l,2})l[g^{t,^2)f^{t,l)^{t,2). 

n=i 12=1 

If the kernel of the interaction potential is the Dirac measure ^{q) = S{q), then the sufficient 
equation for the description of pure state evolution governed by kinetic equation (I150p is the 
Gross-Pitaevskii-type equation 

^ X /* 

t-,p{t, q) = --A,i/j{t, q) + J dq'dq'Xt, q, q; q' , q")i/j{t, q")r{t, q)i/j{t, q'), 

where the coupling ratio b(t,q,q;q',q") is the kernel of the scattering amplitude operator 
Yl^j^=iQi{—t,ii)bi{{l,2})Yl^^^-^^Qi{t,i2). Observing that in the kinetic (macroscopic) scale of 
the variation of variables the groups of operators ([7]) of finitely many particles depend on mi- 
croscopic time variable e~^t, where £ > is a scale parameter, the dimensionless marginal 
functionals of the state are represented in the form: Fs[e~^t,Y \ Fi(t)). Then in the limit 
e — ^ we obtain the Markovian kinetic evolution with the coefficient b(oo, q, q; q', q"). 
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7 Conclusion and outlook 

The concept of cumulants fl6Tl) of the groups of operators forms the basis of the nonperturbative 
solution expansions for hierarchies of quantum evolution equations, namely the dual quantum 
BBGKY hierarchy for marginal observables ( l57l) . the quantum BBGKY hierarchy for marginal 
density operators ( 173|) . the von Neumann hierarchy ( l33l) for correlation operators and the 
nonlinear quantum BBGKY hierarchy for marginal correlation operators ( l92i) . and as well 
as it underlies of the kinetic evolution fllOSp . The nonperturbative hierarchy solutions are 
represented in the form of an expansion over the particle clusters which evolution is governed 
by the corresponding order cumulant (1611) of the groups of operators ([7]) of the Heisenberg 
equations (jl]) or the groups of operators (IT^ of the von Neumann equations ( ITT]) . 

We emphasize that intensional Banach spaces for the description of states of infinite-particle 
systems, which are suitable for the description of the kinetic evolution or equilibrium states, 
are different from the exploit spaces [1], Thus, marginal density operators or correlation 
operators from the space -C^(J-'-h) describe finitely many quantum particles. In order to describe 
the evolution of infinitely many particles we have to construct solutions for initial data from 
more general Banach spaces than the space of sequences of trace class operators. For example, 
it can be the space of sequences of bounded translation invariant operators which contains the 
marginal density operators of equilibrium states In that case every term of the solution 
expansion (177|) of the quantum BBGKY hierarchy or the nonlinear quantum BBGKY hierarchy 
fllUll) and correspondingly the generalized quantum kinetic equation flllOp as well as mean-value 
functional f l54p in case of the dual quantum BBGKY hierarchy fl63p contains the divergent 
traces, which can be renormalized due to the cumulant structure of the solution expansions. 

The origin of the microscopic description of non-equilibrium correlations of Bose and Fermi 
many-particle systems was considered in section 3. For the correlation operators that give 
an alternative description of the quantum state evolution of Bose and Fermi many-particle 
systems, the von Neumann hierarchy of nonlinear evolution equations flS^ was introduced. In 
particular, it was established that in case of absence of correlations in the system at initial 
time, the correlations generated by the dynamics of a system (1381) are completely determined 
by cumulants (13 9 P of the groups of operators (ITS!) of the von Neumann equations. 

The links of constructed solution of the von Neumann hierarchy both with the solution of 
the quantum BBGKY hierarchy (I72p and with the nonlinear quantum BBGKY hierarchy for 
marginal correlation operators ( 15^ were discussed. The cumulant structure of solution f lUp of 
the von Neumann hierarchy fl53l) induces the cumulant structure of solution expansions both 
the initial-value problem of the nonlinear quantum BBGKY hierarchy for marginal correlation 
operators (llOip and the quantum BBGKY hierarchy for marginal density operators (1771) . Thus, 
the dynamics of infinite-particle systems is generated by the dynamics of correlations. We note 
that along with the definition within the framework of the non-equilibrium grand canonical 
ensemble the marginal density operators can be defined within the framework of dynamics of 
correlations (172|) that allows to give the rigorous meaning to the states for more general classes 
of operators than trace-class operators. 

We remark that the rigorous results on the evolution equations in functional derivatives for 
generating functionals of states and observables of classical many-particle systems, namely the 
BBGKY hierarchy and the dual BBGKY hierarchy in functional derivatives respectively, are 
presented in |15] . 
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We developed also an approach to a description of the evolution of states by means of the 
quantum kinetic equations. It was demonstrated that in fact, if initial data is completely 
defined by a one-particle density operator, then all possible states of infinite-particle systems 
at arbitrary moment of time can be described within the framework of a one-particle density 
operator without any approximations. In other words, for mentioned states the evolution of 
states governed by the quantum BBGKY hierarchy (1731) can be completely described by the 
generalized quantum kinetic equation fll03p . It should be emphasized that the kinetic evolution 
is an inherent property of infinite-particle systems. In spite of the fact that in terms of a one- 
particle marginal density operator from the space of trace-class operators can be described a 
system with the finite average number of particles, the generalized quantum kinetic equation 
has been derived on the basis of the formalism of nonequilibrium grand canonical ensemble 
since its framework is adopted to the description of infinite-particle systems in suitable Banach 
spaces as well. 

We note that constructed marginal functionals of the state fll06p or flll5p characterize the 
correlations of states of quantum many-particle systems. Owing to that from macroscopic point 
of view the evolution of many-particle states with the minimal dispersion is the Markovian ki- 
netic evolution, then from microscopic point of view such evolution is characterized by marginal 
correlation functionals (11151) which equal to zero. 

An approach to the description of kinetic evolution of quantum many-particle systems in 
terms of the evolution of marginal observables is also developed. In the mean field limit the 
evolution of marginal observables is governed by the dual quantum Vlasov hierarchy (11251) . One 
of the advantages of such approach as well as developed approach of the generalized quantum 
kinetic equation (11031) is the possibility to construct the kinetic equations in scaling limits 
in case of the presence of correlations at initial time (I150p . for instance, correlations which 
characterize the condensate states of particles [3]. 
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